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Phase  change  material  slurries  (PCS)  can  be  employed  in  the  refrigeration  and  air  conditioning  systems 
as  both  secondary  refrigerant  and  cold  energy  storage  media  simultaneously,  which  benefits  not  only 
the  system  efficiency  improvement  by  the  high  cold  carry  capacity  but  also  the  reduction  of  the 
environment-negative-impact  gas  emission.  This  paper  reviews  the  previous  researches  and  develop¬ 
ments  on  the  applications  of  PCS,  including  ice  slurry,  microencapsulated  phase  change  material  slurry 
(MPCS)  as  well  as  clathrate  hydrate  slurry  (CHS),  to  the  secondary  loop  refrigeration  and  air 
conditioning  systems.  The  paper  mainly  focuses  on  the  generation  approaches  and  storage  strategies 
of  these  slurries  which  are  considered  as  the  most  important  issues  relevant  to  the  application  of  PCS, 
while  the  flow  and  heat  transfer  characteristics  of  PCS  in  both  tubes  and  heat  exchangers  are  also 
summarized.  Moreover,  several  application  cases  of  PCS  and  the  corresponding  operation  performances 
are  presented. 
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1.  Introduction 

The  application  of  the  secondary  refrigerant  to  the  refrigeration 
and  air  conditioning  systems  can  reduce  the  charging  amount  of 
the  environment-negative-impact  refrigerant,  like  CFCs  and  HCFCs, 
etc.,  consequently  the  emission  of  these  refrigerants  to  the  atmo¬ 
sphere  is  reduced.  Meanwhile,  the  cold  energy  storage  is  often 
employed  in  the  refrigeration  and  air  conditioning  systems  to  shift 
the  peak-load  to  the  off-peak  time,  and  thus  many  merits  can  be 
achieved:  reducing  the  system  scale,  minimizing  the  partial  load 
compressor  operation  to  improve  its  efficiency,  improving  the 
system  COP  by  the  extended  night  time  operation  at  lower 
environment  temperature,  reducing  the  financial  cost  by  producing 
and  storing  the  cold  energy  during  night  time  when  the  electricity 
is  cheaper  and  releasing  the  cold  energy  during  day  time  when  the 
electricity  is  more  expensive.  The  inspiring  aspect  is  that  both  the 
roles  of  the  secondary  refrigerant  and  cold  energy  storage  medium 
can  be  undertaken  by  phase  change  material  slurries  (PCS)  simul¬ 
taneously,  such  as  ice  slurry,  microencapsulated  phase  change 
material  slurry  (MPCS)  and  clathrate  hydrate  slurry  (CHS).  The 
latent  heat  of  the  fusion  of  the  solid  particles  containing  in  PCS 
makes  it  possible  to  preserve  cold  energy  with  high  cold  storage 
density,  therefore  the  volumes  or  sizes  of  the  utilized  storage  tank, 
piping,  pump  and  the  refrigerator  are  relatively  small;  the  easily 
pumpable  performance  of  PCS  through  tubes  as  well  as  heat 
exchangers  without  clogging  provides  the  qualification  for  PCS 
being  a  good  secondary  refrigerant  with  high  cold  carry  capacity, 
then  the  pumping  power  is  significantly  reduced  since  the  required 
flow  rate  is  low.  Therefore,  PCS  was  considered  promising  as 
“functional  thermal  fluids”  [1]. 

The  utilization  of  ice  slurry  has  been  widely  investigated  over 
the  world  from  the  fundamental  thermal  properties  to  the 
practical  applications  [2-4].  Ice  slurry  is  a  homogeneous  mixture 
of  very  small  ice  crystals,  usually  0.1-2  mm  in  characteristic  size, 
and  the  base  liquid  can  be  either  pure  water  or  binary  solution  of 
water  and  freezing  point  depressing  agent,  e.g.  ethanol,  glycol  or 
salt,  etc.  The  latent  heat  of  ice  crystals  is  about  333.55  kj  kg 
which  results  in  a  high  cold  carry  capacity.  Compared  to  single¬ 
phase  water  or  glycol  solution,  ice  slurry  with  30  wt%  ice  fraction 
offers  an  80%  reduction  of  flow  rate  and  therefore  the  pumping 
power  while  maintaining  the  same  cooling  load  [5].  Due  to  its 
excellent  cooling  and  transporting  features,  ice  slurry  can  be 
employed  in  many  fields,  such  as  the  building  cooling,  super¬ 
market  display,  food  industry,  brewery,  mine  cooling,  fish  pro¬ 
cessing,  as  well  as  medical  cooling.  However,  the  practical 
application  is  still  limited  by  the  following  disadvantages:  the 
evaporation  temperature  is  low  for  air  conditioning  using;  heavy 
mechanical  power  is  required  to  generate  small  and  smooth  ice 
crystals  and  to  maintain  the  ice  slurry  in  homogeneous  state;  the 
pressure  drop  and  pumping  power  of  ice  slurry  with  high  ice 
fraction  are  very  high.  Therefore,  the  utilization  of  ice  slurry  is  not 


common  except  in  Japan,  and  only  few  projects  exist  in  USA  and 
Europe  up  to  date. 

MPCS  is  a  kind  of  solid-liquid  suspension  consisting  of  the  carrier 
liquid  and  small  particles  of  phase  change  material  (PCM)  with  a 
thin  shell.  The  diameter  of  the  microcapsules  is  in  the  range  of 
1-1000  pm,  the  core  PCM  can  be  solid,  liquid,  or  even  mixture. 
Paraffin  waxes  are  often  used  as  the  core  materials  due  to  its 
extended  melting  temperature  range.  The  shell  material  is  usually 
natural  or  synthesized  polymer,  such  as  polyester  and  polyethylene, 
which  has  high  strength  and  flexibility.  The  carrier  liquid  is  usually 
water.  Compared  to  non-encapsulated  PCM  slurry,  MPCS  owns 
several  advantages:  the  agglomeration  of  PCM  particles  and  the  in 
turn  clogging  in  piping  are  avoided;  the  PCM  does  not  contact  with 
environment  directly,  which  avoids  disadvantageous  changes  of 
PCM  composition  and  property;  the  microencapsulation  offers  a 
protective  structure  element  which  improves  the  mechanical  stabi¬ 
lity,  thus  permits  the  frequent  changes  in  volume  of  the  core 
material  as  phase  change  occurs;  large  heat  transfer  area  is  achieved 
by  the  surface  of  the  microencapsulation.  So  far,  however,  the 
investigations  about  MPCS  are  mostly  on  its  fabrication  methods 
and  thermal  properties  [6-8]  and  the  flow  and  heat  transfer 
characteristics  through  tubes  [9-15],  Due  to  the  high  cost  of  MPCS 
preparation,  only  few  laboratorial  applications  of  MPCS  to  chilled 
ceiling  panels  [16,17]  and  an  industrial  cooling  project  at  Narita 
Airport  in  Japan  [18]  are  reported  in  the  open  literature. 

Clathrate  hydrate  can  be  formed  in  some  aqueous  solution  of 
tetraalkylammonium  salts  with  simple  anions  (such  as  halides, 
sulfate,  formate,  etc.)  [19]  at  the  atmospheric  condition,  in  which 
the  water  molecules  form  structure  lattices  as  host  and  molecules 
of  another  substance  fill  the  lattices  as  guest;  also,  gases,  like  C02, 
N2,  CFCs  and  so  on,  can  form  hydrate  in  water  under  the  higher 
pressure  and  low  temperature  conditions;  however,  such  type  of 
gas  hydrate  slurry  would  not  be  discussed  in  the  present  paper 
because  its  application  require  more  rigorous  conditions. 
Researchers  in  Japan  [20]  applied  CHS  to  the  air  conditioning 
system  as  cold  storage  and  transportation  medium,  and  found 
that  tetrabutylammonium  bromide  (TBAB)  CHS  was  quite  suita¬ 
ble  to  the  air  conditioning  application  due  to  its  following  merits: 
TBAB  CHS  can  be  easily  generated  at  atmospheric  pressure  and 
room  temperature  by  cooling  down  the  aqueous  solution;  TBAB 
CHS  has  a  freezing  point  between  0  and  12  °C  depending  on  TBAB 
concentration  in  the  aqueous  solution,  which  fits  quite  well  with 
the  necessary  coolant  temperature  in  the  air  conditioning  system; 
the  cooling  capacity  of  TBAB  CHS  is  about  2-4  times  as  large  as 
that  of  chilled  water  within  the  temperature  range  of  5-12  °C;  the 
fluidity  of  TBAB  CHS  is  also  good.  Researchers  in  JFE  Engineering 
Corporation  also  successfully  applied  TBAB  CHS  to  a  4000  m2 
office  building  [21]  for  cooling  and  the  operational  results  showed 
47%  reduction  of  the  pumping  power.  Except  the  application 
projects  conducted  by  JFE  Engineering  Corporation  [21,22];  however, 
other  studies  of  TBAB  CHS  still  focused  on  the  fundamental 
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Nomenclature 

Cp  specific  heat  (kj  kg-1  K-1) 

D  diameter  (m) 

dp  particle  size  (m) 

/  friction  factor  (dimensionless) 

fF  Fanning  friction  factor  (dimensionless) 

Gz  Graetz  number  (dimensionless),  Gz=RePrDhlx 

h  heat  transfer  coefficient  (W  m~2  I<_1) 

AH  enthalpy  change  (kj  kg-1) 

H  height  (m) 

K'  fluid  consistency  coefficient  (dimensionless) 

L  length  (m) 

n  flow  behavior  index  (dimensionless) 

Nu  Nusselt  number  (dimensionless) 

Pr  Prandtl  number  (dimensionless) 

AP  pressure  drop  (Pa) 

qw  wall  heat  flux  (W) 

Re  Reynolds  number  (dimensionless) 

ReMR  modified  Reynolds  number  (dimensionless) 

Ste  Stefan  number  (dimensionless) 

T  temperature  (°C) 

AT  temperature  difference  (°C) 

v  flow  velocity  (m  s~]) 

xs  particle  fraction  (wt%) 


Axs  particle  fraction  difference  (wt%) 

x  distance  from  inlet  (m) 

Greeks  symbols 

/ (  viscosity  (mPa  s) 

X  thermal  conductivity  (W  m_1  K-1) 

cp,  <P  volume  fraction  (dimensionless) 

Subscripts 

a  average  value 

B  Bingham 

BT  Bingham  parameters  in  Tomita  equation 

cf  carrier  fluid 

/  fluid 

h  hydraulic 

i  inner 

is  ice  slurry 

local  local  value 

m  mean  value 

overall  overall  value 

s  slurry 

w  wall 


properties  [23],  the  flow  and  heat  transfer  characteristics  of  TBAB 
CHS  flowing  through  tubes  [24-31]  as  well  as  in  heat  exchangers 
[32,33],  In  addition,  other  CHS  materials,  such  as  trimethylo- 
lethane  (TME)  hydrate  slurry  [34-37]  and  tetrabutylammonium 
fluoride  (TBAF)  hydrate  slurry  [38],  have  also  been  studied  as  the 
thermal  energy  storage  medium. 

This  paper  presents  a  review  of  the  current  main  applications 
of  PCS,  including  the  aforementioned  ice  slurry,  MPCS  and  TBAB 
CHS,  to  the  refrigeration  and  air  conditioning  systems  as  both  the 
cold  storage  medium  and  secondary  refrigerant.  The  generation 
methods  and  storage  strategies  which  are  very  important  to  PCS 
application  are  summarized,  and  the  flow  and  heat  transfer 
characteristics  of  PCS  in  tubes  and  heat  exchangers  are  also 
discussed.  Some  application  cases  of  the  practical  projects  or 
laboratorial  studies  are  displayed,  and  the  corresponding  opera¬ 
tion  results  are  also  discussed. 


2.  Generation  approaches  and  storage  strategies 

2.1.  Generation  approaches  of  ice  slurry  and  TBAB  CHS 

The  generation  processes  of  ice  crystals  in  water  or  aqueous 
solution  and  that  of  TBAB  clathrate  hydrate  are  similar  in  that  the 
crystals  initiate  and  grow  up  in  the  solution.  The  generation 
process  generally  include  the  following  steps:  supercooling  of  the 
base  liquid,  nucleation  and  growth  of  the  crystals.  The  base  liquid 
is  first  cooled  down  to  the  temperature  lower  than  its  phase 
change  temperature  but  without  phase  change,  the  state  of  which 
is  named  as  the  supercooled  state.  Then  the  nucleation  occurs, 
which  is  usually  heterogeneous  since  the  alien  dust  particles  as 
well  as  the  rough  container  wall  provide  nuclei  for  the  nucleation. 
And  then  it  is  followed  by  the  growth  of  crystals.  After  these  three 
steps,  size  and  shape  of  the  crystals  are  affected  by  the  processes  of 
attrition,  agglomeration  and  ripening  [39].  If  the  generated  slurry  is 
agitated  or  collides  with  other  materials,  attrition  will  happen  due 
to  shear  stress  on  the  crystals,  and  large  crystals  will  break  off 


during  the  attrition.  Agglomeration  may  occur  during  the  slurry 
generation  or  storage  without  stirring,  and  small  crystals  collide 
and  stick  to  each  other  to  form  the  large  size  crystals.  If  the  slurry  is 
stored  for  long  time,  small  crystals  may  dissolve  slowly  and  then 
the  size  distribution  changes,  which  is  called  Ostwald  ripening. 

Usually,  small  and  globular  crystals  are  expected  for  easy 
pumping  through  tubes  as  well  as  heat  exchangers,  therefore 
the  attrition  process  is  commonly  utilized  for  crystal  size  reduc¬ 
tion  and  edge  smoothing.  On  the  other  side,  it  was  reported  that 
the  ice  crystals  produced  in  pure  water  were  dendritic  and  had 
very  rough  surface,  which  resulted  in  the  poor  fluidity  and  easy 
agglomeration  and  clogging  even  though  the  crystal  size  was 
small  [4[.  The  improvement  on  this  dendritic  structure  can  be 
achieved  by  smoothing  the  crystals  in  chemical  way  through 
adding  the  freezing  point  depressing  agents  to  water,  such  as 
ethanol,  glycol,  propylene  glycol  and  sodium  chloride  and  so  on. 
Hayashi  and  Kasza  [40]  produced  ice  slurry  in  glycol  solution,  and 
found  that  even  very  small  amount  of  glycol  could  improve  the 
fluidity  of  the  generated  ice  slurry  significantly,  and  the  higher 
concentration  the  better  fluidity. 

The  generation  approaches  are  the  pivotal  technology  for  the 
applications  of  both  ice  slurry  and  CHS,  since  not  only  the  crystal 
properties,  i.e.,  size,  shape  and  roughness,  are  determined  by  the 
generation  process,  but  also  the  system  energy  efficiency  and 
reliability  mainly  lie  in  the  generator  performance.  The  genera¬ 
tion  of  ice  slurry  has  been  widely  investigated  in  the  alternative 
ways,  which  can  be  mainly  categorized  to  the  methods  using 
moving  parts  and  the  methods  without  moving  parts.  Therefore, 
the  following  contexts  are  mainly  focused  on  the  research  results 
of  ice  slurry,  and  these  approaches  can  also  be  applied  to  CHS 
because  of  the  similar  generation  mechanism  as  aforementioned. 


2.1.1.  Scraped  surface  generator 

Scraped  surface  ice  slurry  generator  has  been  widely  investi¬ 
gated  and  developed  over  the  world,  and  it  has  been  applied  in 
large-scale  systems  [5],  The  scraped  surface  generator  commonly 


5024 


P.  Zhang,  Z.W.  Ma  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  5021-5058 


utilizes  a  shell-and-tube  heat  exchanger  in  which  the  evaporating 
refrigerant  flows  through  the  shell  side  while  the  crystals  are 
generated  in  the  tube  side.  Driven  by  a  motor,  a  rotatable  scraper 
is  mounted  on  the  center  of  the  inner  tube  to  prevent  the  crystals 
from  sticking  to  the  cold  heat  transfer  surface,  as  shown  in 
Fig.  1(a)  [39],  Stamatiou  [41]  and  Stamatiou  and  Kawaji  [42] 
developed  a  compact  heat  exchanger  with  scraper  inside  the 
channel  to  generate  ice  slurry,  as  shown  in  Fig.  1(b).  The  scraper 
moves  up  and  down  to  eliminate  the  ice  crystals  on  the  tube  wall. 
The  authors  obtained  30%  power  saving  by  using  compact  heat 
exchanger  instead  of  tubular  one,  since  four  times  lower  agitation 
speed  was  required  to  maintain  the  ice-free  heat  transfer  surface. 

The  blades  of  the  scraper  may  be  plastic  or  brush  type  so  that 
they  can  have  good  contact  with  the  tube  wall  to  detach  the 
crystal  layer,  since  it  is  believed  that  the  nucleation  took  place  on 
the  coldest  area  inside  the  tube,  i.e.,  the  inner  tube  wall,  and  then 
the  crystals  layer  grows  on  the  wall  forming  a  thermal  resistance 
layer  which  reduces  the  heat  transfer  coefficient  between  the 
refrigerant  and  fluid  inside  the  tube.  However,  other  researchers 
thought  that  the  crystals  were  formed  in  the  bulk  of  the  super¬ 
cooled  liquid  but  had  a  tendency  to  adhere  to  the  cold  wall  by  the 
cohesion  force,  and  this  was  why  the  blades  were  needed.  More¬ 
over,  the  scraper  movement  also  contributes  to  the  great  enlarge¬ 
ment  of  the  fluid  turbulent  degree,  therefore  the  heat  transfer 
performance  of  the  heat  exchanger  is  promoted.  Rodriguez 
Pascual  et  al.  [43]  investigated  the  influences  of  the  plastic  blade 
geometries  on  crystal  layer  removal  efficiency  and  the  slurry 
mixing  effect  by  visualization  and  simulation.  The  authors  found 
that  the  mixing  of  scraped  ice  crystals  with  the  bulk  solution  by 


Fig.  1.  Schematic  diagram  of  the  scraped  surface  ice  slurry  generators,  (a)  shell- 
and-tube  type  [39];  (b)  compact  heat  exchanger  type  [41]  and  (c)  auger  type 
scraper  [47  J. 


the  streamlined  blade  was  inadequate,  while  the  vertical  blade 
promoted  the  mixing  but  accumulated  the  crystals  in  front  of  it. 

Alternatively,  Lakhdar  et  al.  [44]  employed  stainless  steel  blades 
rather  than  the  plastic  ones.  The  authors  permitted  a  small  gap 
between  the  blade  and  the  tube  wall,  and  found  the  permanent  ice 
layer  only  existed  when  the  gap  was  large  enough,  about  1  mm 
according  to  their  study.  This  fact  may  indicate  that  the  ice  crystals 
nucleate  inside  the  solution.  Helical-blade  auger  can  also  be 
utilized  as  the  scraper  for  ice  slurry  generation  [45-47],  as  shown 
in  Fig.  1(c).  In  addition  to  the  effect  of  scraping  the  heat  transfer 
surface,  the  rotating  auger  also  provides  a  negative  pumping 
pressure  to  suck  the  water  (or  solution)  flowing  through  the 
generator,  hence  it  can  provide  the  function  of  the  pump. 

The  generation  of  ice  slurry  using  scraped  surface  is  considered 
as  very  reliable,  even  in  the  case  of  little  freezing  point  depressing 
agent  being  added  to  water.  During  the  generation,  the  increase  of 
ice  fraction  is  usually  implemented  by  adjusting  the  liquid  flow 
rate.  Using  a  Sunwell  scraped  surface  ice  slurry  generator,  Wang 
and  Kusumoto  [5]  produced  0-35  wt%  ice  slurry  with  about 
100  pm  ice  crystals,  and  the  used  base  liquid  is  an  aqueous 
solution  with  6-10%  ethylene  glycol.  However,  drawbacks  exist 
in  this  generator:  the  tight  contact  of  the  plastic  scraper  blade  on 
the  tube  wall  causes  strong  friction,  and  then  the  blades  need  to 
be  replaced  after  a  certain  time  interval  of  utilization;  since  the 
generator  employs  additional  mechanical  components,  the  initial 
investment  and  the  energy  consumption  are  high. 

The  scraped  surface  generator  is  also  suitable  to  TBAB  CHS 
generation.  Douzet  et  al.  [48]  successfully  applied  a  scraped 
surface  heat  exchanger  (about  200  L,  supplied  by  HeatCraft  Ltd., 
France)  to  TBAB  CHS  air  conditioning  system.  Daitoku  and  Utaka 
[49]  investigated  the  scrape  effect  of  different  scraper  shapes  on 
the  detachment  of  TBAB  clathrate  hydrate  from  heat  transfer 
surface.  The  authors  found  that  type  B  TBAB  clathrate  hydrate 
cohered  to  the  surface  more  tightly  and  needed  stronger  force  to 
detach,  and  a  relief  angle  of  the  scraper  was  required  to  avoid  the 
accumulation  of  crystals  on  the  clearance  between  the  scraper 
and  the  heat  transfer  surface  when  detaching  type  B  hydrate. 


2.1.2.  Orbital  rod  ice  slurry  generator 

Another  type  of  ice  slurry  generator  using  moving  parts  is  the 
orbital  rod  generator,  as  shown  in  Fig.  2  [50].  This  generator  is 
usually  a  vertical  shell-and-tube  heat  exchanger  with  evaporating 
refrigerant  flowing  from  bottom  to  top  in  the  shell  side,  while  inside 
the  shell  is  a  tube  bundle  consisted  many  thin  tubes.  As  shown  in 
the  figure,  each  of  tube  has  a  rotatable  whip  rod  inside,  which  does 
not  rotate  in  the  tube  center  but  around  the  tube  circumference.  The 
rotating  rods  undertake  several  roles:  pushing  part  of  the  upper 
income  liquid  flowing  around  the  tube  circumference  to  form  a 
falling  film;  enhancing  the  film  coefficient  and  the  flow  turbulence 
and  therefore  the  heat  transfer  coefficient  between  the  film  and  the 
refrigerant;  preventing  the  generated  ice  crystals  from  sticking  to 
the  heat  transfer  surface,  just  like  the  function  of  scraper  in  the 
scraped  surface  ice  slurry  generator.  The  rotation  of  the  whip  rod  is 
much  faster  (e.g.,  850  rpm)  than  that  of  the  scraper  (e.g.,  450  rpm), 
however,  it  is  believed  that  the  falling  film  acts  as  the  lubricant 
which  minimizes  the  drive  power  of  the  motor. 

Orbital  rod  ice  slurry  generator  is  also  widely  used,  e.g.,  Paul 
Mueller  Company  [51]  developed  ice  slurry  units  using  this 
generator  which  has  the  standard  production  sizes  from  3  to 
520  t  (about  10.55-1828.68  kW).  They  utilized  the  generator  in 
many  fields,  including  refrigeration,  HVAC  and  process  cooling  and 
so  on.  Stamatiou  et  al.  [39]  compared  the  performances  of  the 
orbital  rod  and  the  scraped  surface  ice  slurry  generators.  It  is 
believed  that  both  the  two  kinds  of  generators  are  stable,  and  the 
generated  ice  crystals  are  homogeneous  due  to  influence  of  the 
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Fig.  2.  Schematic  diagram  of  the  orbital  rod  ice  slurry  generator  [50], 


rapid  rotated  scraper  or  rods.  The  scraped  surface  ice  slurry 
generator  has  a  slight  higher  cooling  capacity,  however,  the  drive 
power  for  the  scraper  rotation  is  also  larger;  the  overall  heat 
transfer  coefficient  of  the  orbital  rod  ice  slurry  generator  is  higher 
due  to  the  falling  film  heat  exchange.  The  whip  rod  is  metallic  and 
has  no  direct  contact  with  the  tube,  therefore  it  survives  longer 
time  than  that  of  the  plastic  scraper  blade  which  attaches  to  the 
tube  tightly  as  aforementioned.  However,  the  drive  mechanism  of 
whip  rods  is  complicated  and  there  are  many  associated  problems. 


2.1.3.  Fluidized  bed  generator 

The  fluidized  bed  heat  exchanger  was  proposed  in  1970s  for 
sea  water  desalination  [52].  Small  solid  particles  are  inserted  into 
the  tubes  in  a  shell-and-tube  heat  exchanger  to  detach  the 
adherence  on  the  tube  wall.  The  insert  particles  are  usually 
cylindrically  shaped  stainless  steel  made  by  wire,  or  globular 
glass  balls.  As  the  fluid  flowing  upward  in  the  tube,  the  solid 
particles  are  fluidized  and  continuously  impact  on  the  heat 
transfer  wall,  therefore  the  deposits  on  the  wall  are  broken  off 
and  fall  into  the  fluid,  as  shown  in  Fig.  3  [53].  Meanwhile,  the  heat 
transfer  coefficient  is  also  improved  by  the  following  three 
mechanisms  [53]:  the  laminar  boundary  layer  that  acts  as  heat 
transfer  resistance  layer  is  broken  up  by  the  fluidized  particles; 
the  solid  particles  transport  the  heat  from  the  heat  transfer  wall 
to  the  fluid  by  conduction  as  the  particles  contacting  the  wall  and 
then  the  fluid;  the  deposit  layer  on  the  wall  is  broken  off. 
Therefore,  the  heat  transfer  coefficient  is  determined  by  the  flow 
velocity  of  the  fluid  and  also  the  bed  voidage  (the  ratio  of  the 
space  that  not  occupied  by  the  solid  particles  to  the  entire  space). 

Klaren  and  van  der  Meer  [54]  first  applied  the  fluidized  bed 
heat  exchanger  to  produce  ice  slurry  in  1991,  then  Meewisse  and 
Infante  Ferreira  [55,56]  and  Pronk  et  al.  [53,57-59]  conducted  a 
lot  of  works  on  producing  ice  slurry  by  the  fluidized  bed. 
Schematic  diagrams  of  the  typical  fluidized  bed  ice  slurry  gen¬ 
erators  are  shown  in  Fig.  4.  The  generator  shown  in  Fig.  4(a)  is  the 
stationary  fluidized  bed,  in  which  the  insert  solid  particles  are 
kept  inside  the  tube  during  the  operation.  Shown  in  Fig.  4(b)  is  the 
so-called  circulating  fluidized  bed,  where  the  solid  particles  will 


Fig.  3.  Schematic  diagram  of  the  effect  of  fluidized  solid  particles  [53], 

flow  out  of  the  tube  together  with  the  generated  ice  slurry,  and 
then  are  separated  from  the  slurry  and  return  to  the  bottom  via  a 
downcomer.  Measurement  results  by  Meewisse  and  Infante 
Ferreira  [56]  showed  that  the  heat  transfer  coefficients  of  the 
fluidized  bed  during  ice  slurry  generation  was  in  the  range  of 
2500-4000  W  m~2  K  \  which  was  about  2-4  times  as  large  as 
that  of  single-phase  flow  and  was  not  influenced  by  the  ice  slurry 
generation.  They  also  proposed  the  empirical  heat  transfer  corre¬ 
lations  which  predicted  the  heat  transfer  accurately.  Compared  to 
the  stationary  mode,  the  use  of  the  circulating  fluidized  bed  is 
more  flexible  since  the  solid  particle  can  flow  at  a  high  velocity 
range,  and  the  advantages  include  the  improvement  of  the  ice 
crystals  breaking  off  from  the  tube  wall  by  increasing  flow 
velocity  and  also  the  heat  transfer  coefficient  is  increased  by 
1-25%  [58], 
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Fig.  4.  Schematic  diagram  of  the  fluidized  bed  ice  slurry  generator,  (a)  stationary  type  and  (b)  circulating  type  [57], 


The  shell-and-tube  design  of  the  fluidized  bed  is  simple  and 
results  in  an  inexpensive  and  low  operation  cost  for  ice  genera¬ 
tion.  Due  to  the  violent  impacts  of  solid  particles  on  the  tube  wall, 
the  generated  ice  particles  can  be  homogeneous,  which  is  usually 
in  the  shape  of  cylindrical  or  ellipsoidal  with  the  size  of  about 
300  pm  [57],  The  main  disadvantage  of  the  fluidized  bed  ice 
slurry  generator  is  the  existence  of  a  maximum  allowable  tem¬ 
perature  difference  between  the  fluid  inside  the  tube  and  the 
refrigerant  (or  the  tube  wall),  below  which  ice  slurry  can  be 
generated  stably  with  nearly  constant  heat  transfer  coefficient. 
The  higher  temperature  difference  enlarges  the  amount  of  ice 
crystals  sticking  on  the  tube  wall,  which  are  hard  to  be  removed 
efficiently  by  the  impact  of  solid  particles.  The  gradually  thicken¬ 
ing  crystal  layer  deteriorates  the  heat  transfer  between  the  slurry 
and  the  refrigerant,  thereafter  ice  slurry  generation  becomes 
unstable.  Meewisse  and  Infante  Ferreira  [55]  measured  this 
temperature  difference  limit  in  a  circulating  fluidized  bed,  as 
shown  in  Fig.  5.  It  was  found  that  the  temperature  difference  limit 
varied  linearly  with  the  freezing  temperature  of  the  aqueous 
solution.  Lower  freezing  temperature  (or  higher  additive  concen¬ 
tration)  caused  a  higher  temperature  difference  limit.  The  other 
disadvantage  of  this  ice  slurry  generator  is  the  requirement  of 
higher  fluid  velocity,  so  that  the  momentums  of  the  solid  particles 
are  sufficient  to  impact  the  wall  for  ice  crystal  layer  detachment. 


2.1.4.  Supercooling  method 

Ice  slurry  generation  using  supercooling  concept  employs  the 
following  steps:  water  or  aqueous  solution  is  cooled  in  an 
evaporator  to  a  few  degrees  below  its  freezing  temperature 
without  crystallization;  the  supercooled  liquid  is  forced  to  flow 
through  a  special  device,  i.e.,  supercooling  releaser,  to  release  the 
supercooling  state;  then  ice  crystals  are  generated  in  the  liquid 
and  the  liquid  turns  to  solid-liquid  ice  slurry.  This  method  is 
considered  as  simple  and  having  low  operation  cost  since  there  is 
no  other  complex  apparatus.  However,  the  evaporator  for  super¬ 
cooling  and  the  supercooling  releaser  must  be  carefully  designed 
to  achieve  stable  ice  slurry  generation. 

The  heat  exchanger  for  supercooling  must  ensure  that  there  is 
no  crystallization  happened  inside,  otherwise  the  slurry  genera¬ 
tion  may  stop  due  to  the  blockage  by  the  accumulatively  gener¬ 
ated  ice  crystals.  Bedecarrats  et  al.  [60]  regarded  that  the 
stochastic  character  of  crystallization  resulted  in  a  hard  control 
of  ice  generation  within  heat  exchanger  and  it  was  necessary  to 


Fig.  5.  Range  of  stable  ice  slurry  production  using  circulating  fluidized  bed. 
(a)  range  of  NaCl  and  (b)  comparison  of  the  operating  ranges  of  NaCl  and  ethylene 
glycol  (EG).  (A Tmax  is  the  temperature  difference  between  the  slurry  and  the 
refrigerant,  the  bed  voidage  is  0.8)  [55], 


control  the  evaporation  temperature  accurately.  They  studied  the 
influences  of  the  supercooling  degree,  refrigerant  temperature 
and  fluid  flow  rate  on  the  ice  generation.  By  developing  models 
for  predicting  the  fluid  temperature  variation,  they  tried  to  avoid 
the  ice  generation  inside  supercooling  heat  exchanger. 
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Several  requirements  for  the  supercooling  heat  exchanger  are 
given  as  follows: 

(1)  The  heat  transfer  surface  should  be  smooth,  since  the  rough 
surface  enhances  the  crystallization.  Through  the  observations  of 
crystallization  within  tubes  with  different  roughness  (between 
0.63  and  13.3  pm),  Faucheux  et  al.  [61]  found  strong  influence  of 
surface  roughness  on  supercooling  degree  and  concluded  that 
larger  roughness  caused  lower  supercooling  degree. 

(2)  The  flow  channel  for  the  water  or  solution  should  be  regular, 
since  the  sudden  change  of  the  channel,  e.g.  bend,  contraction 
or  expansion,  results  in  strong  collision  between  liquid  and 
solid  surface,  which  accelerates  the  nucleation. 

(3)  The  heat  exchange  inside  the  evaporator  should  be  uniform 
everywhere  of  the  heat  exchanger;  otherwise  crystallization 
may  occur  locally  where  the  temperature  of  the  fluid  is  low. 

The  supercooling  releaser  is  also  important  to  the  effective 
crystal  generation,  and  several  releasing  approaches  are  investi¬ 
gated  and  utilized  as  follows: 

(1)  The  collision  of  the  supercooled  fluid  with  a  solid  surface,  or 
free  liquid  surface  or  liquid  bulk  in  the  tank,  or  with  itself,  as 
shown  in  Fig.  6  [62-65].  The  collision  provides  the  momen¬ 
tum  increase  to  release  the  supercooling  state.  The  probability 
of  freezing  mainly  depends  on  the  collision  force  [66].  More¬ 
over,  the  system  requires  special  device  to  avoid  ice  crystals 
agglomeration  near  the  impact  place. 

(2)  Ultrasonic  vibration :  Researchers  [67,68]  reported  that  the  ultra¬ 
sonic  vibration  can  strongly  promote  the  transition  of  super¬ 
cooled  water  to  ice  slurry,  and  this  method  is  reliable  for  ice 
generation  but  it  is  only  limited  to  laboratorial  investigation. 

(3)  Partial  low  temperature:  Partial  liquid  being  cooled  to  the 
temperature  much  lower  than  its  freezing  temperature  can  be 


helpful  to  the  nucleation.  This  process  can  be  achieved  by  a 
small  electronic  refrigerator. 

(4)  Electric  charge:  Though  scientists  have  tried  to  investigate  the 
mechanism  of  the  nucleation  by  electric  charge,  the  basic 
mechanism  is  still  controversial  [69].  However,  some  experi¬ 
mental  results  have  been  identified  to  shown  the  influential 
factors.  Okawa  and  Saito  [66]  disclosed  supercooling  released  by 
the  utilization  of  electric  charge  with  the  voltage  less  than  DC 
100  V,  and  found  that  the  influential  factors  were  the  voltage,  the 
size  of  electrode  and  the  distance  between  the  two  electrodes. 

Okawa  and  Saito  [66]  also  treated  the  reciprocally  rubbing  of 
glass  against  glass  as  one  of  the  effective  supercooling  release 
approach,  which  was  depended  upon  the  frequency  and  ampli¬ 
tude  of  the  rubbing.  On  the  other  hand,  the  authors  believed  the 
convection,  vibration,  and  non-contacting  shock  had  no  effect  on 
supercooling  release. 

Tanino  and  Kozawa  [70]  and  Kozawa  et  al.  [71]  applied  a  pipe 
type  releaser  for  large-scale  ice  slurry  generation,  which  generated 
2.5  wt%  ice  slurry  for  one  circulation,  and  they  finally  obtained 
35  wt%  ice  slurry.  While  Ise  et  al.  [72]  used  this  kind  of  releaser  to 
generate  about  23  wt%  ice  slurry  for  a  building  usage. 

Nearly  the  same  approach  can  be  employed  for  CHS  generation, 
e.g.,  JFE  Engineering  Corporation  [73]  successfully  applied  TBAB  CHS 
air  conditioning  to  several  large-scale  buildings  [21,22,74].  Ogishi 
et  al.  [75]  employed  a  tube  type  heat  exchanger  (such  as  shell-and- 
tube  heat  exchanger)  as  the  supercooling  evaporator  for  TBAB  CHS 
generation.  The  heat  exchanger  also  took  the  role  of  supercooling 
releaser,  therefore  there  was  at  least  one  of  the  tubes  had  irregularity 
for  supercooling  release.  Hydrates  would  be  generated  inside  the  heat 
exchanger.  The  authors  reported  the  deposition  of  hydrate  crystals  on 
the  tube  inner  wall  and  treated  them  to  be  useful  nuclei  for  stable 
hydrate  generation,  while  the  other  hydrates  were  removed  by  the 
shear  force  of  the  flowing  aqueous  solution  or  slurry.  The  thickness  of 


Fig.  6.  Collision  modes  for  supercooling  release,  (a)  supercooled  liquid  and  solid  surface  [63];  (b)  supercooled  liquid  and  liquid  surface  in  the  tank  [64];  (c)  supercooled 
liquid  and  bulk  liquid  in  the  tank  [65]  and  (d)  supercooled  liquid  and  itself  [62], 
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the  deposition  layer  was  further  controlled  by  adjusting  the  sluny 
flow  rate  and  regulating  the  refrigerant  temperature.  Moreover,  Song 
et  al.  [76]  validated  the  effect  of  ultrasonic  vibration  on  TBAB  CHS 
generation  in  25  wt%  TBAB  solution,  and  the  nucleation  was  observed 
to  be  generated  in  a  very  short  time  under  the  ultrasonic  vibration 
and  the  crystallization  process  also  finished  in  several  seconds.  The 
authors  found  that  the  crystals  may  be  cracked  under  the  force  of 
ultrasonic  vibration,  and  the  obtained  TBAB  hydrate  crystals  were 
fine,  smooth  and  uniform  without  agglomeration. 


2.1.5.  Direct  contact  heat  exchanger 

Using  the  direct  contact  heat  transfer  between  water  or  aqueous 
solution  with  immiscible  cool  refrigerant,  including  primary  and 
other  secondary  refrigerants,  ice  slurry  can  be  generated  directly 
inside  the  storage  tank  [77-83].  Fig.  7(a)  shows  the  utilization  of 
primary  refrigerant.  During  the  evaporation  of  the  two-phase  primary 
refrigerant,  dispersed  ice  particles  are  generated  while  the  refrigerant 
can  be  easily  separated  since  it  turns  to  vapor.  However,  the  heat 
exchanger  needs  to  work  under  high  pressure.  Fig.  7(b)  shows  the 
application  of  secondary  refrigerant  which  has  lower  freezing  tem¬ 
perature  and  usually  larger  density  than  water,  the  refrigerant  sinks 
to  bottom  of  the  tank  after  heat  exchange  and  then  returns  to  the 
pump  and  chiller  for  another  cycle. 
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Fig.  7.  Schematic  diagrams  of  the  direct  contact  heat  exchanger  for  ice  slurry 
generation,  (a)  using  primary  refrigerant  [78]  and  (b)  using  secondary  refrigerant  183], 


Within  a  storage  tank  of  0.3  m  in  diameter  and  1.1  m  in  height, 
Kiatsiriroat  et  al.  [78]  utilized  the  primary  refrigerant,  including 
R22  and  R12  in  the  temperature  range  of  -8  to  -15°C,  as  the 
coolant  for  direct  contact  heat  exchange  with  water.  The  flow  rates 
of  the  two  refrigerants  were  0.02-0.08  kg  s-1  (about  0.28- 
1.13  kg  m-2  s-1)  while  water  was  stagnant  and  kept  inside  the 
tank.  They  found  that  the  volumetric  heat  transfer  coefficients  in 
the  cases  of  R12  and  R22  were  around  2-16kWm-3K-1  and 
2-52  kW  m  3  K-1,  respectively.  The  promoted  heat  transfer  per¬ 
formance  by  R22  was  attributed  to  its  higher  water  solubility. 
Thitipatanapong  et  al.  [79]  investigated  ice  generation  using  R12 
and  R134a  in  direct  contact  heat  exchange  while  both  the  refrig¬ 
erant  and  water  were  forced  flowing  into  the  heat  exchanger.  It 
was  noticed  that  higher  water  flow  rate  resulted  in  higher  volu¬ 
metric  heat  transfer  coefficient.  The  volumetric  heat  transfer 
coefficients  obtained  by  the  authors  was  about  91.62  kW  m-3  K  1 
and  209.7  kW  m-3  K_1  in  the  cases  of  R12  and  R134a  with  water 
flow  rates  of  about  16.35  kg  m-2  s-1  and  16.20  kg  m  2  s_1, 
respectively,  which  was  higher  than  those  obtained  by  Kiatsiriroat 
et  al.  [78].  Due  to  the  saturation  pressure  of  refrigerant,  the  ice 
slurry  forming  temperature  increased  to  4.0  °C  and  8.8  °C  in  the 
case  of  R12  and  R134a,  respectively,  therefore  the  storage  tem¬ 
perature  at  about  4.5  °C  and  10.0  °C  can  be  achieved,  which  was 
considered  suitable  to  air  conditioning  system  application.  Thong- 
wik  et  al.  [80]  injected  low  temperature  C02  (- 15  to  -60  °C)  to 
water  at  a  flow  rate  of  0.003-0.017  kgs-1,  and  the  corresponding 
volumetric  heat  transfer  coefficient  was  ranged  from  1  to 
25  kW  m-3  K  After  the  direct  contact  with  the  primary  refrig¬ 
erant,  the  base  liquid  is  inevitably  contaminated  by  the  compressor 
oil,  therefore  the  authors  added  the  compressor  oil  and  surfactant 
Tween  60  (polyoxyethylene  sorbitan  monostearate,  which  was 
used  to  merge  water  and  oil)  to  water  to  investigate  its  influence 
on  sluny  generation.  It  was  found  lower  oil  concentration 
increased  the  risk  of  blockage  around  injector  since  the  non¬ 
merging  water  formed  plain  ice  and  blocked  the  injector,  and  the 
optimum  composition  of  water/oil/surfactant  was  100/6/1. 

Researchers  in  Singapore  [81-83]  proposed  a  special  coolant, 
named  Fluroinert  (FC-84),  as  the  suitable  secondary  refrigerant  used 
in  direct  contact  ice  slurry  generation.  FC-84  is  heavier  than  water, 
the  density  of  which  is  about  1730  kg  m-3,  while  the  freezing 
temperature  is  -  95  °C.  As  shown  in  Fig.  7(b),  FC-84  was  sprayed  by 
a  nozzle  at  the  top  of  the  storage  tank.  The  coolant  was  at  about 
-10  to  -15°C  while  the  flow  rate  was  8-12  L  min-1.  Through 
modeling  and  calculation,  the  overall  heat  transfer  coefficient  of 
the  heat  exchange  between  water  and  FC-84  was  in  the  range  of 
3.0-6.5  kW  m  2  l<  '.  The  generated  ice  crystals  were  in  the  shape 
of  thin  and  long  flake  with  the  length  of  about  1  -4  mm,  which  was 
quite  large  compared  to  those  obtained  by  other  methods. 

Compared  to  other  ice  slurry  generators,  the  utilization  of 
direct  heat  exchange  avoids  the  sticking  of  ice  crystals  on  the  heat 
transfer  surface  and  strongly  enhances  the  heat  transfer  perfor¬ 
mance.  Therefore  the  COP  of  the  refrigerator  increases,  for 
example,  about  3.4-3. 6,  higher  than  that  of  the  conventional 
system  [77].  The  disadvantages  of  this  method  are  as  follows:  the 
nozzle  used  in  the  system  must  be  well  designed,  otherwise 
blockage  may  occur;  small  amount  of  primary  refrigerant  trapped 
in  the  ice  crystals  will  be  released  to  the  environment. 


2.3.6.  Vacuum  ice  slurry  generator 

Ice  slurry  can  be  generated  in  a  deep  vacuum  freezer.  Water 
inside  the  freezer  is  kept  at  deep  vacuum  close  to  the  triple-point 
of  water,  hence  part  of  the  water  becomes  vapor,  which  absorbs 
the  latent  heat  of  vaporization  from  the  left  water,  and  it  is  then 
frozen  and  converted  to  ice  particles  or  ice-water  mixture. 
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Fig.  8.  Vacuum  ice  slurry  generator  used  in  thermal  energy  storage  system  [85], 

IDE  Technologies  Ltd.  [84-86]  developed  a  vacuum  ice  slurry 
generator,  and  the  schematic  diagram  is  shown  in  Fig.  8.  A  high 
speed  centrifugal  compressor  driven  by  a  motor  was  utilized  to 
generate  deep  vacuum  condition.  The  generated  ice  slurry  inside 
the  freezer,  containing  0.5-1. 0  mm  ice  crystals,  was  then  pumped 
to  the  storage  tank,  and  the  slurry  fraction  could  reach  50%.  To 
generate  ice  slurry  continuously,  water  vapor  was  compressed 
and  then  forced  flowing  through  a  condenser,  where  the  vapor 
was  cooled  by  chilled  water  (e.g.,  5  °C)  and  converted  to  liquid 
phase.  The  cooling  capacity  of  this  generator  can  achieve 
3500  kW.  The  generator  has  been  applied  to  district  cooling  and 
heating,  thermal  energy  storage,  mine  cooling  and  snow  making, 
and  hundreds  units  have  been  installed  worldwide  [86].  However, 
due  to  large  specific  volume  of  water  vapor  under  deep  vacuum 
condition,  the  size  of  the  generator  is  quite  large.  On  the  other 
side,  Bapat  and  Kulkarni  [87]  employed  a  lithium  bromide  vapor 
absorption  machine  to  generate  the  vacuum  condition  for  ice 
slurry  generation,  which  is  compact  in  size  and  less  electric 
energy  consumption. 

Researchers  in  Korea  [88,89]  sprayed  water  droplets  with 
initial  temperature  at  20  °C  and  50  pm  in  size  by  nozzles  in  a 
vacuum  chamber  to  generate  ice  slurry.  The  pressure  in  the 
chamber  was  kept  lower  than  the  pressure  of  triple  point,  and 
spherical  ice  crystals  with  size  below  300  pm  were  then  formed. 
To  avoid  the  generation  of  ice  crystals  around  nozzles,  the  authors 
added  7%  ethylene  glycol  to  water.  Different  from  the  vacuum 
generator  shown  in  Fig.  8,  the  formed  vapor  was  then  removed 
from  the  chamber.  The  feeding  of  water  was  maintained  at  low 
rate  to  obtain  water  droplets,  therefore  the  cooling  capacity  was 
only  0.23  kW  (theoretical  value),  which  could  be  improved  by 
increasing  the  nozzle  number. 

2.2.  Storage  strategy  of  ice  slurry  and  TBAB  CHS 

The  storage  of  ice  slurry  and  TBAB  CHS  in  tanks  is  mainly 
conducted  in  two  modes,  homogeneous  storage  and  heteroge¬ 
neous  storage  [90],  as  shown  in  Fig.  9. 

As  shown  in  Fig.  9(a),  the  homogeneous  storage  required 
agitation  to  maintain  the  slurry  fluidizing  inside  the  storage  tank 
and  to  avoid  the  agglomeration  of  particles.  The  agitation  is 
usually  provided  by  propeller,  pumps,  paddles  and  so  on.  Since 
the  absence  of  agglomeration  in  homogeneous  slurry  storage,  it  is 
very  easy  to  extract  the  slurry  with  very  fine  crystals  and  the 
pumping  is  also  quite  reliable.  It  was  reported  [91]  that  the 
maximum  ice  fraction  in  homogeneous  storage  could  be  up  to 
30-35  wt%,  beyond  which  the  power  consumed  by  the  agitator 
increased  significantly,  e.g.,  the  power  consumed  by  Danish 


Fig.  9.  Storage  strategies  of  ice  slurry,  (a)  homogeneous  storage  and 
(b)  heterogeneous  storage  [90], 


Technological  Institute  to  agitate  52  wt%  ice  slurry  was  10  times 
larger  than  that  of  30  wt%. 

Since  the  agitation  is  significant  energy  consuming,  Meili  et  al. 
[92]  reported  that  the  mixing  of  the  slurry  could  be  stopped  for 
about  1 2  h,  and  then  homogeneous  ice  slurry  could  be  still 
obtained  after  restarting  the  agitation.  Moreover,  adding  freezing 
point  depressing  agents  to  water  can  smooth  the  generated  ice 
particles  and  prevent  the  agglomeration.  However,  Grandum  and 
Nakagomi  [93]  and  Inaba  et  al.  [94]  believed  the  additions  of 
freezing  point  depressing  agents  resulted  in  a  low  refrigerator 
COP  because  of  the  depressed  freezing  temperature.  Therefore, 
the  authors  utilized  other  additives,  for  example,  antifreeze 
proteins  (AFPs),  antifreeze  glycoproteins  (AFPs),  silane  coupling 
agents,  surfactants  and  so  on  and  considered  them  as  suitable 
since  the  low  additive  concentration  (less  than  1  wt%)  led  to 
relative  high  freezing  temperature. 

Fig.  9(b)  shows  the  method  of  heterogeneous  storage.  Due  to 
the  different  density  of  the  solid  crystals  from  that  of  the  carrier 
liquid,  such  as  ice  is  lighter  than  water  and  TBAB  hydrate  is 
heavier  than  the  solution,  crystals  in  heterogeneous  storage  tank 
tend  to  accumulate  at  the  top  or  deposit  at  the  bottom  of  the  tank. 
Compared  to  homogeneous  storage,  the  ice  fraction  in  hetero¬ 
geneous  storage  tank  can  reach  about  50-60  wt%  [91  ],  and  the 
necessary  volume  of  the  storage  tank  is  then  reduced.  The  difficult 
thing  in  this  storage  mode  is  that  the  extraction  of  ice  from  the 
tank.  Eglof  et  al.  [90]  considered  using  a  screw  pump  to  extract 
the  high  fraction  ice  slurry  from  the  top  of  the  tank  and  then 
mixed  it  with  the  low  fraction  ice  slurry  pumped  from  the  bottom 
region  through  a  special  valve.  This  proposal  has  not  been 
experimentally  validated  yet  and  needs  further  investigation. 
Douzet  et  al.  [48]  employed  this  approach  for  TBAB  CHS  storage, 
they  pumped  the  high  fraction  TBAB  CHS  (about  40-50%  in 
volume)  from  the  bottom  of  the  tank  using  a  lobes  pump  and 
mixed  it  with  the  liquid  sucked  through  the  top  of  the  tank  using 
another  pump,  then  they  obtained  TBAB  CHS  with  about  20% 
volume  fraction  for  air  conditioning  application.  Actually,  pure 
ice  particles  can  be  obtained  by  a  harvester  [91],  though  the 
shaved  ice  particles  (e.g.,  1-2  mm)  were  larger  than  the  original 
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generated  ones  (e.g.,  0.1  mm),  the  particles  are  mixed  with  the 
cold  water  extracted  from  the  bottom  of  the  tank  to  form  ice 
slurries  with  different  fractions. 

2.3.  Fabrication  of  MFCS 

MPCS  requires  the  technology  that  coating  the  PCM  with 
suitable  shell  material  to  envelop  the  core  PCM.  Gibbs  et  al.  [95] 
reported  several  forms  of  the  microcapsules,  as  shown  in  Fig.  10. 
Generally  three  steps  are  required  for  microcapsules  preparation: 
forming  the  shell  around  the  PCM,  ensuring  there  is  no  leakage, 
and  ensuring  the  undesired  materials  are  not  included  in  the 
capsule.  Though  the  microencapsulation  technology  was  discov¬ 
ered  and  developed  in  the  1940s  and  1950s  [96]  and  then  many 
approaches  were  developed  for  food  industrial  [97],  only  coacer- 
vation,  polymerization,  spray  drying  were  found  being  employed 
in  PCM  encapsulation.  These  processes  are  simply  introduced  in 
the  following  contexts,  the  detailed  information  about  the  fabri¬ 
cation  and  the  further  properties  of  the  formed  microencapsu¬ 
lated  PCMs  can  refer  to  the  review  articles  [7,8]. 

2.3 A.  Coacervation 

There  are  two  types  of  coacervations:  simple  coacervation  and 
complex  coacervation.  The  complex  one  is  commonly  used  for 
PCM  microencapsulation.  First,  polymer  in  water  as  the  coating 
material,  e.g.  gelatine  in  water,  is  prepared.  Then  the  core 
material  is  dispersed  to  the  gelatine  solution  at  a  temperature 
range  of  40-60  C,  at  which  the  coating  material  is  liquid.  The 
simple  coacervation  employs  only  one  kind  of  colloid,  while 
another  colloid,  such  as  gum  Arabic,  is  added  to  the  solution  for 
complex  coacervation  at  this  step.  The  step  is  followed  by 
coacervation.  Simple  coacervation  is  induced  by  change  of  condi¬ 
tion,  such  as  the  temperature,  the  addition  of  nonsolvent  or 
microions,  which  causes  the  dehydration  of  macromolecules. 
Complex  coacervation  is  usually  induced  by  increasing  the  poly¬ 
mer  concentration  or  reducing  the  solution  pH  (usually  4-4.5). 
The  system  is  then  cooled  to  room  temperature.  Finally,  the 
microcapsules  are  hardened  and  isolated.  Hardening  is  conducted 
by  cross-linking  agent,  such  as  formaldehyde,  the  increasing  of  pH 
to  9-11  (by  NaOH  solution,  e.g.)  and  temperature  to  50  °C  can 
hasten  this  process.  The  temperature  should  be  reduced  and 


maintained  at  about  5-10  °C  for  2-4  h,  then  the  capsules  are 
filtered  and  dried  in  a  vacuum  oven. 

Using  the  complex  coacervation,  Hawlader  et  al.  [98]  obtained 
around  1  pm  spherical  and  smooth  microcapsules  of  paraffin 
wax.  The  latent  heat  was  about  234.05Jg~\  224.31  Jg-1  and 
196. 38  J  g_1  when  the  core-to-coating  ratio  was  2:1,  1:1  and  1:2, 
respectively.  The  encapsulation  efficiency  was  in  the  range  of 
80-95%,  and  it  decreased  as  increasing  the  ratio  of  core-to-coating 
since  incomplete  coating  on  core  material  might  occur  when  the 
coating  material  was  deficient.  Ozonur  et  al.  [99]  applied  the 
same  method  to  encapsulate  the  natural  coco  fatty  mixture  for 
thermal  energy  storage  (melting  temperature  of  the  PCM  was  in 
the  range  of  22-24  °C).  The  size  of  the  obtained  microcapsules 
varied  from  1  pm  to  1  mm  depending  on  the  stirring  speed  and 
time.  After  encapsulation,  there  was  an  increment  of  about  7  °C 
on  the  melting  temperature. 

2.3.2.  Polymerization 

Polymerization  are  used  widely  for  fabricating  microencapsu¬ 
lated  PCM,  and  there  are  several  different  type  polymerization 
approaches.  Yang  et  al.  [100]  and  Brown  et  al.  [101]  employed 
the  so-called  “in-situ”  polymerization  to  encapsulate  tetradecane 
and  dicyclopentadiene,  respectively.  The  former  study  utilized 
several  types  of  polymers,  e.g.,  polyvinyl  acetate  (PVAc),  polystyrene 
(PS),  polymethylmethacrylate  (PMMA),  polyethyl  methacrylate 
(PEMA).  The  authors  mixed  the  PCM  and  polymer  with  surfac¬ 
tants,  initiator  and  deionized  water  in  a  flask  accompanying  with 
agitation,  and  then  heated  the  mixture  to  and  maintained  it  at  a 
certain  temperature.  The  mixture  was  then  cooled  down  naturally 
by  continuous  stirring  for  4-5  h.  It  was  found  that  the  coating 
material  strongly  influenced  the  thermal  properties  of  the  formed 
MPCM.  After  encapsulation,  the  melting  temperature  decreased 
from  6.7  °C  (pure  tetradecane)  to  5.97  °C  as  the  coating  material 
was  PMMA  and  to  5.68  °C  as  the  coating  material  was  PEMA, 
while  the  freezing  temperature  increased  from  1.98  °C  to  2.95  C 
and  to  3.19  °C,  respectively.  The  obtained  microcapsules  were 
5-30  pm  in  diameter  and  66.26  kj  kg-1  in  melting  heat  with  40% 
tetradecane.  The  latter  study  described  their  detailed  preparation 
process  in  Fig.  11.  Using  this  process,  10-1000  pm  microcapsules 
were  obtained  by  adjusting  the  agitation  rate  between  200  and 
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Fig.  10.  Several  forms  of  the  microcapsules  [95]. 


Fig.  11.  In-situ  polymerization  for  microencapsulation  [101]. 
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2000  rpm,  and  the  wall  thickness  was  in  the  range  of  160- 
220  nm.  It  was  noticed  that  the  surface  roughness  reduced  as 
the  agitation  rate  increased.  Fill  content  of  the  obtained  micro¬ 
capsules  was  beneficially  up  to  83-92  wt%.  Using  in-situ  poly¬ 
merization  and  the  assistance  of  ultrasonic,  Fang  et  al.  [6]  coated 
n-octadecane  with  polystyrene  as  shell  and  obtained  very  fine 
nanoencapsulated  PCM,  the  particles  were  regular  spheres  and 
diameter  was  only  50-200  nm.  The  phase  change  temperature 
was  found  very  close  to  that  of  pure  n-octadecane,  while  the 
latent  heat  dropped  nearly  50%  due  to  the  existence  of  the  coating 
material. 

Alkan  et  al.  [102]  prepared  the  microcapsules  using  emulsion 
polymerization.  They  used  docosane  as  core  and  PMMA  as  shell. 
Docosane  was  first  mixed  with  deionized  water  and  Triton  X-100 
(surfactant),  and  then  was  heated  to  the  temperature  higher  than 
the  melting  temperature  of  docosane  under  nitrogen  atmosphere. 
This  process  was  followed  by  the  addition  of  methylmethacrylate 
and  allyl  methyl  acrylate  and  FeS04  ■  7H20  solution,  of  which  the 
allyl  methyl  acrylate  was  used  as  cross-linking  agent  to  prevent 
the  formation  of  secondary  particles.  The  mixture  was  then 
agitated  violently  for  30  min.  Na2S207  solution  and  tertbutylhy- 
droperoxide  solution  were  added  and  the  mixture  was  heated  to 
90  C.  Thereafter  the  polymerization  started.  After  casting  water, 
the  precipitate  was  dried  under  vacuum  at  40  °C  for  24  h.  The 
average  diameter  of  the  obtained  PMMA  microcapsules  was 
0.16  pm.  The  influences  of  microencapsulation  on  the  melting 
and  freezing  temperature  were  slight,  and  only  0.6  C  decrease  of 
the  melting  temperature  and  0.2  °C  increase  of  the  freezing 
temperature  were  observed  after  the  microencapsulation.  The 
melting  and  freezing  heat  were  measured,  which  were 
54.6  kjkg-1  and  48.7  kjkg-1,  respectively.  Therefore  the  slurry 
contained  28  wt%  docosane. 

Cho  et  al.  [103]  and  Siddhan  et  al.  [104]  applied  almost  the 
same  interfacial  polymerization  to  encapsulate  the  PCM  (octade- 
cane)  using  TD1  (toluene-2,  4-diisocyanate)  and  DETA  (diethyle- 
netriamine)  as  coating  material.  Surfactant  NP-10  [poly( ethylene 
glycol)  nonylphenyl  ether]  was  treated  as  emulsifier  by  the 
former  study  while  the  latter  study  used  NP-9.5.  The  NP-10  (or 
NP-9.5)  aqueous  solution  and  organic  solution  of  octadecane, 
cyclohexane  and  TDI  were  prepared  and  mixed  with  each  other 
accompanying  with  violent  agitation  to  form  an  oil-in-water 
emulsion.  Then  the  DETA  aqueous  solution  was  added  to  the 
emulsion  drop  by  drop  to  start  the  interfacial  polymerization 
which  happened  between  TDI  and  DETA  at  the  oil-water  inter¬ 
face.  The  rotation  speed  of  the  agitator  was  slowed  down  and  the 
mixture  was  heated  to  60  DC.  This  state  was  maintained  for 
45  min  (the  latter  study  added  more  DETA  solution  for  further 
diffusion  and  polymerization).  Finally,  the  microcapsules  were 
filtered  and  washed  and  dried.  Both  studies  obtained  small 
capsules,  0.1-1  pm  in  the  former  study  and  less  than  14  pm  in 
the  latter  one.  By  additionally  adding  DETA  solution,  Siddhan 
et  al.  [104]  obtained  higher  core  fill  content  (about  40-68  wt%) 
than  that  obtained  by  Cho  et  al.  (26.0-46.4  wt%).  The  encapsula¬ 
tion  efficiency  was  75.4-87.2%  in  the  former  study,  while  the 
latter  study  noticed  that  the  efficiency  was  largely  depended  on 
the  core-to-monomer  ratio  (varied  from  about  50%  to  90%  as  the 
ratio  increased  from  about  2  to  6.5)  and  PCM-to-cyclohexane 
ratio  (varied  from  about  87%  to  97%  as  the  ratio  increased  from 
about  1  to  12). 

Ai  et  al.  [105]  used  suspension  polymerization  for  n-hexade- 
cane  microencapsulation,  while  polystyrene  was  used  as  shell. 
Considering  casein  in  solution  and  dispersed  form  as  an  excellent 
emulsifier,  the  authors  utilized  casein  instead  of  the  necessary 
surfactants  as  emulsifier  in  other  polymerization  processes  to 
emulsify  the  mixture  of  core  material,  monomers  and  initiators. 
The  usage  of  casein  provided  extra  effect  of  protecting  the  oil 


droplets  against  recoalescence  before  the  microcapsules  forma¬ 
tion.  Thus,  the  application  of  casein  in  polymerization  process 
could  simplify  the  process  since  no  more  additives  were  needed. 
The  size  of  the  obtained  microcapsules  was  3-15  pm.  The  latent 
heat  was  in  the  range  of  76.01-124.7  kj  kg~'  depending  on  the 
polymerization  time,  which  was  much  lower  than  that  of  the  pure 
material,  236kjkg~1.  Sanchez  et  al.  [106]  also  considered  this 
method  as  simple,  cheap  and  technically  feasible.  They  coated 
several  kinds  of  phase  change  materials,  including  paraffin  wax, 
tetradecane,  nonadecane,  etc.,  by  a  coating  material  of  styrene. 
The  other  necessary  materials  included  benzoyl  peroxide  as 
initiator,  polyvinylpyrrolodone  as  stabilizer  and  water.  After  the 
violent  agitation,  the  mixture  was  heated  to  110  °C  at  one  atmo¬ 
sphere  for  6  h,  then  it  was  washed,  filtrated  and  dried.  The 
capsules  size  ranged  between  20  pm  and  700  pm  with  the 
average  value  of  237  pm,  which  was  also  dependent  on  different 
core  materials.  The  core  fill  content  for  paraffin  wax,  tetradecane 
and  nonadecane  was  20.56  wt%,  39.10  wt%  and  49.36  wt%, 
respectively. 

2.3.3.  Spray  drying 

Spray  drying  is  a  technology  that  can  generate  powder 
instantaneously  by  atomizing  the  liquid  production  in  a  hot  air 
or  nitrogen  [107],  which  is  considered  as  economic  for  micro¬ 
encapsulation.  The  spray  drying  process  in  microencapsulation 
involves  the  following  basic  steps:  mixing  the  core  solution  and 
the  coating  material  solution;  heating  and  agitating  the  mixture 
to  obtain  homogeneous  mixture;  if  needed,  adding  the  emulsifier 
to  the  mixture,  which  depends  on  the  properties  of  the  coating 
material;  ensuring  the  formed  emulsion  are  stable  over  a  certain 
period;  atomizing  the  emulsion  in  the  drying  chamber  where  the 
solvent  (usually  water)  is  evaporated;  cooling  and  dehydrating 
the  obtained  particles.  Hawlader  et  al.  [98]  employed  a  spray 
dryer  to  spray  dry  the  emulsions  of  paraffin  wax  and  the  coating 
material  with  a  feeding  rate  of  20mLmin~1  by  a  centrifugal 
atomizer  operated  at  25,000  rpm.  The  authors  compared  the 
microcapsules  produced  by  spray  drying  with  that  by  complex 
coacervation,  and  it  was  found  the  spray  drying  made  much 
smaller  particles  (around  0.1  pm  compared  to  around  1  pm), 
however,  the  latent  heat  was  also  smaller  than  that  produced 
by  the  complex  coacervation  (about  221.5  kjkg-1  compared  to 
234  kj  kg-1,  when  the  core-to-coating  ratio  was  2:1). 

3.  Flow  and  heat  transfer  characteristics 

Flow  and  heat  transfer  characteristics  of  PCS  in  tubes  and  heat 
exchangers  are  significant  to  its  application.  The  rheology,  flow 
and  heat  transfer  behaviors  of  ice  slurry  have  been  previously 
reviewed  by  Ayel  et  al.  [108]  in  2003,  and  Zhang  et  al.  [7]  also 
presented  the  flow  and  heat  transfer  characteristics  of  MPCS  and 
CHS  in  2010.  However,  there  are  many  new  researches  about  the 
PCS  that  have  been  conducted  thereafter.  Therefore  we  emphasize 
the  present  review  on  the  recent  research  details  of  the  flow 
pattern,  flow  frictions  and  heat  transfer  coefficients  of  PCS  in 
tubes  and  heat  exchangers  and  also  the  impact  of  particle  fraction 
and  tube  diameter  on  mentioned  behaviors. 

3.1.  Flow  pattern  of  ice  slurries 

The  usual  flow  patterns  of  solid-liquid  two-phase  flow  are 
homogeneous  flow,  heterogeneous  flow,  moving  bed  and  station¬ 
ary  bed.  The  pressure  drop  increases  with  flow  velocity  at  the 
homogeneous  and  heterogeneous  flows  but  decreases  at  moving 
bed  and  stationary  bed  flows.  All  these  flow  patterns  have  ever 
been  observed  in  ice  slurry  flow,  however,  Kitanovski  et  al.  [109] 
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pointed  that  the  flow  pattern  classification  by  simple  observation 
showed  quite  large  uncertainty. 

Using  the  Ultrasonic  Doppler  Profiler  (UVP),  Vuarnoz  et  al.  [110] 
obtained  both  the  vertical  and  horizontal  velocity  profiles  of 
7.5-7. 6%  ice  slurry  flowing  through  a  20  mm  x  20  mm  rectangular 
channel,  as  shown  in  Fig.  12(a).  With  the  decrease  of  flow  rate, 
both  velocity  profiles  became  bulgy  in  the  axes.  The  vertical 
profiles  were  more  sensitive  to  the  flow  rate  reduction  due  to  the 
gravity  effect,  therefore  the  flow  became  slower  at  the  top  and 
faster  at  the  bottom  since  the  accumulation  of  ice  particles  at  the 
top,  as  shown  in  the  figure;  meanwhile,  the  horizontal  profiles 
remain  symmetrical  as  decreasing  flow  rate.  Using  a  hot  film 


anemometer  (HFA)  and  an  on-line  ice  slurry  sampling/calorimetry 
technique,  Stamatiou  and  Kawaji  [111,112]  obtained  the  velocity 
profiles  of  0-5.9  vol%  ice  slurry  and  3. 1-9.0  vol%  ice  slurry  with 
heat  exchange  flowing  through  a  vertical  610  mm  x  310  mm  x 
25.4  mm  channel  as  shown  in  Fig.  12(bl)-(b4).  When  ice  fraction 
was  low,  i.e.  0-0.9  vol%  in  Fig.  12(bl),  the  velocity  profiles  were 
hardly  influenced  by  the  ice  particles,  therefore  the  results  were 
similar  to  that  of  single-phase  fluid;  ice  particles  influence  the 
velocity  profiles  when  the  ice  fraction  is  in  the  range  of  0-5.9  vol%, 
as  shown  in  Fig.  12(b2),  larger  ice  fraction  results  in  blunter 
profiles;  when  heat  flux  was  imposed  to  ice  slurry  flow,  the 
formation  of  ice-free  layer  near  the  heated  wall  led  to  larger  fluid 
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Fig.  12.  Velocity  profiles  of  ice  slurry  flow,  (a)  xs=7.5-7.6%,  rectangular  channel  20  mm  x  20  mm  [110];  (bl)-(b4)  base  liquid  was  6.2  wt%  NaCl  brine,  cfm=0.1-0.2  mm, 
vertical  rectangular  channel  610  mm  x  310  mm  x  25.4  mm  [111,112],  (bl)  adiabatic  flow  at  very  low  ice  fraction;  (b2)  adiabatic  flow  for  various  ice  fractions;  (b3)  diabatic 
flow  for  various  heating  powers;  (b4)  diabatic  flow  for  various  ice  fractions.  The  reference  and  equations  in  (bl)-(b4)  are  those  in  [111,112]. 
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density  in  this  area  than  fluid  in  the  flow  center,  therefore  the  flow 
near  the  heated  wall  decelerated  as  shown  in  Fig.  12(b3)  while  the 
flow  in  the  tube  center  changed  little;  Fig.  12(b4)  shows  that  there 
is  no  clear  influence  of  ice  fraction  on  the  velocity  profiles  of  ice 
slurry  under  heating. 

Ice  fraction  profile  of  ice  slurry  is  related  to  many  factors,  e.g., 
density  of  carrier  liquid,  density  of  ice  particle,  shape  and  size  of 
ice  particle,  tube  diameter,  and  the  flow  velocity  [109].  Kitanovski 
and  Poredos  [113]  modeled  and  numerically  calculated  the 
concentration  profiles  of  10wt%  and  25  wt%  ice  slurry  flow 
though  27.2  mm  tube,  as  shown  in  Fig.  13.  It  was  believed  that 
the  10  wt%  ice  slurry  flow  was  homogeneous  at  14=2.0  m  s_1  and 
was  heterogeneous  at  us=0.5-1.5  m  s_1  and  moving  bed 
appeared  at  vs=0.5  m  s_1,  as  shown  in  Fig.  13(a),  and  the  larger 
ice  faction  resulted  in  earlier  transition  from  homogeneous  flow 
to  heterogeneous  flow  as  comparing  Fig.  13(a)  and  (b). 

With  the  reduction  of  flow  velocity,  flow  of  ordinary  solid- 
liquid  two-phase  slurry  will  transit  from  heterogeneous  flow  to 
moving  bed  at  a  critical  velocity,  which  is  named  as  deposition 
velocity.  On  the  basis  of  velocity  profiles,  Kitanovski  [114] 
evaluated  the  deposition  velocities  of  ice  slurry  with  different 
ice  particle  diameters  as  well  as  different  tube  diameters,  as 
shown  in  Fig.  14(a).  As  the  ice  fraction  increasing  from  5  wt%  to 
25  wt%,  the  deposition  velocity  of  ice  slurry  had  a  maximum  value 
at  xs=15  wt%.  Both  the  increase  of  tube  diameter  and  ice  particle 
diameter  caused  easier  deposition  corresponding  to  the  larger 
deposition  velocity  as  shown  in  the  figure.  Rached  et  al.  [115] 
reported  a  new  method  to  determine  the  deposition  velocity  of 
ice  slurry.  They  measured  the  cross-section  areas  of  the  moving 
bed  as  increasing  the  flow  velocity  gradually,  and  then  extra¬ 
polated  the  deposition  velocity  corresponding  to  the  disappear¬ 
ance  of  the  moving  bed  which  could  not  been  observed  directly. 
Deposition  velocities  of  10-32  wt%  ice  slurry  flowing  through 
50  mm  and  75  mm  tubes  were  obtained,  as  shown  in  Fig.  14(b). 
As  seen  from  the  figure,  the  deposition  velocity  reduced  with  ice 
fraction  but  increases  with  tube  diameter. 

On  the  other  hand,  different  from  the  aforementioned  flow 
patterns,  Hirochi  et  al.  [116]  observed  three  types  of  flow  patterns 
of  ice  slurry  with  chipped  ice  particles.  The  authors  believed  that 
the  behavior  of  ice  slurry  was  different  from  that  of  ordinary 
solid-liquid  two-phase  slurry  due  to  the  cohesive  nature  of  ice 


particles.  The  cohesive  nature  of  ice  particles  dominated  ice  slurry 
flow  pattern  while  the  gravity  dominated  that  without  cohesive 
nature.  Therefore,  the  states  of  moving  bed  and  stationary  bed 
were  not  observed  but  instead  of  cluster  flow  and  column  flow, 
where  the  pressure  drop  also  decreased  as  the  reduction  of  the 
flow  velocity,  as  shown  in  Fig.  15. 
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Fig.  14.  Deposition  velocities  of  ice  slurry  flow  in  tubes,  (a)  Kitanovski  [114]  and 
(b)  Rached  et  al.  [115]. 
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Fig.  13.  Ice  fraction  profiles  of  ice  slurry  flow  (simulated),  base  liquid  was  10%  ethanol/water  mixture,  dm= 1  mm,  D,= 27.2  mm,  (a)  xs=10  wt%  and  (b)  xs=25  wt%  [113]. 
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3.2.  Rheology  of  PCS  and  the  friction  factor  correlations 
3.2 A.  Ice  slurry 

It  was  commonly  considered  that  ice  slurry  with  low  ice 
fraction,  i.e.  6-15%  [108],  behaved  as  Newtonian  fluid.  Therefore 
the  flow  friction  can  be  easily  predicted  by  empirical  correlations 
while  the  apparent  viscosity  was  usually  calculated  by  Thomas 
relation  [108,109] 

/i  =  /t/(l+2.5<p  +  10.05(p2  +  0.00273e16-6'?’)  (1) 


U  (m/s) 

a  b  c 


Fig.  15.  Flow  patterns  and  pressure  drops  of  ice  slurry  with  chipped  ice  particles 
in  tube  [116].  (a)  Dispersed  flow,  02.5  m/s,  (b)  cluster  flow,  01.5  m/s  and 
(c)  column  flow,  01.0  m/s. 


Using  Poiseuille  equation  for  laminar  flow  and  Blasius  correlation 
for  turbulent  flow,  Grozdek  et  al.  [117]  concluded  that  the 
calculation  discrepancies  were  smaller  than  10%  for  both  the 
laminar  flow  of  ice  slurry  with  xs  <  20  wt%  in  9  mm  tube  and  the 
turbulent  flow  with  xs<15wt%  in  9  mm  and  15  mm  tubes. 
Fernandez-Seara  et  al.  [118]  also  treated  the  5-15  wt%  ice  slurry 
as  Newtonian  fluid  in  a  fan-coil  unit  and  predicted  the  flow  friction 
with  relative  error  of  -21.5  to  29  %  in  laminar  region  by  Poiseuille 
equation,  while  with  relative  errors  of  -  6.3  to  9%  and  - 1 7  to  1  %  in 
turbulent  region  by  Petukhov  and  Blasius  equations,  respectively. 

However,  many  other  researchers  considered  ice  slurry  as  non- 
Newtonian  fluid  and  developed  semi-empirical  correlations  based 
on  rheological  properties,  e.g.,  the  Bingham  model,  the  Ostwald 
model  (also  known  as  the  Power-law  model),  the  Casson  model 
and  the  Herschel-Bulkley  model. 

Based  on  the  experimental  flow  curves  of  0-30  wt%  ice  slurries 
in  10-20  mm  tubes  as  shown  in  Fig.  16(a),  Niezgoda-Zelasko  and 
Zalewski  [119]  determined  ice  slurry  as  Bingham  fluid,  which 
indicated  that  the  ice  slurries  would  behave  as  Newtonian  fluid 
after  conquering  the  yield  stress.  The  corresponding  yield  stress 
and  apparent  viscosity  for  ice  slurries  with  different  ice  fractions 
were  measured,  then  the  friction  of  laminar  flow  was  estimated  by 
Buckingham-Reiner  equation,  i.e.,  Eq.  (2),  with  excellent  agree¬ 
ment  and  the  friction  factor  of  turbulent  flow  was  calculated  by 
Tomita  equation,  i.e.,  Eq.  (3),  with  the  uncertainty  better  than  10% 


1  _ f p  He  He4 

ReB~  16  61 ?e2  +  3 fFRe8B 
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where  ReBT  was  calculated  by  the  equation  given  by  Tomita  [119]. 

Different  from  Niezgoda-Zelasko  and  Zalewski  [119],  Kumano 
et  al.  [120]  treated  ice  slurry  as  pseudo-plastic  fluid  based  on  their 
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Fig.  16.  Rheological  feature  determination  of  ice  slurry,  (a)  shear  stress  and  shear  rate  of  ice  slurry  flowing  through  16  mm  tube  [1 19];  (b)  flow  curve  of  ice  slurry  flowing 
through  10.2  mm  tube  [120]  and  (c)  apparent  viscosity  of  ice  slurry  measured  by  rheometer  [121]. 
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measured  flow  curves  of  6-24 wt%  ice  slurries  in  4.3-10.2  mm 
tubes,  as  shown  in  Fig.  16(b).  The  authors  utilized  the  Power-law 
model  to  illustrate  the  flow  characteristics  of  ice  slurry,  and 
obtained  the  corresponding  parameters:  flow  behavior  index  n 
and  fluid  consistency  coefficient  K' ,  then  the  modified  Reynolds 
number,  ReMR,  was  used  to  describe  the  flow  of  power-law  ice 
slurries,  as  given  in  the  following: 


ReMR  = 


Dvp 

nr 


(  4n 


\3n  +  l 


Dnv2~np 

K'S"-1 


(4) 


It  was  found  that  the  ice  slurry  deviated  from  Newtonian  fluid 
more  with  higher  ice  fraction  and  in  larger  tube. 

Stokes  et  al.  [121]  measured  the  apparent  viscosities  of 
0-29  vol%  ice  slurries  by  a  rheometer,  results  are  shown  in 
Fig.  16(c).  As  shown  in  the  figure,  there  was  a  large  zero-shear 
viscosity  at  low  applied  stress,  after  which  the  apparent  viscosity 
dropped  drastically  as  the  increasing  shear  stress.  Therefore,  the 
authors  suggested  that  the  Herschel-Bulkley  model  should  be 
applied  to  ice  slurries. 

Grozdek  et  al.  [117]  tried  several  types  of  models  and  friction 
correlations  on  the  flow  of  0-30  wt%  ice  slurry  in  9-25  mm  tubes 
and  revealed  that:  Buckingham-Reiner  equation  for  laminar  flow 
predicted  the  experimental  results  with  the  relative  error  within 
+  15%  for  only  60%  of  all  results,  and  the  inclusion  of  stratification 
effect  would  extend  this  error  range  to  82%  of  all  experimental 
data;  Dodge-Metzner  correlation  for  Power-law  fluid,  i.e.,  Eq.  (5) 
and  Tomita  equation  for  Bingham  fluid  showed  good  agreement 
with  the  experimental  results  in  turbulent  flow 
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However,  Reynolds  numbers  utilized  in  the  study  were  not 
corresponding  to  each  rheological  model,  but  just  used  Reynolds 
number  defined  as  Newtonian  fluid  (apparent  viscosity  was 
calculated  by  Thomas  relation). 

Niezgoda-Zelasko  and  Zelasko  [122]  proposed  a  general 
method  to  estimate  the  flow  friction  of  ice  slurry  flowing  through 
passages  with  arbitrary  cross-section.  The  calculation  process  was 
similar  to  that  of  Power-law  method,  however,  the  geometrical 
parameters  of  the  flow  passage  was  added  to  the  expressions  of 
flow  behavior  index  and  fluid  consistency  coefficient.  ReMR  was 
then  obtained,  using  which  Poiseuille  equation  for  laminar  flow 
and  Dodge-Metzner  correlation  as  well  as  Blasius  equation  for 
turbulent  flow  could  be  employed  by  non-Newtonian  ice  slurry 
flow.  The  comparison  between  the  calculated  and  measured 
friction  factors  of  ice  slurry  flow  through  circular  tube  (16  mm 
in  diameter)  and  rectangular  channel  (3  mm  x  38.5  mm  and 
78  mm  x  26.5  mm)  showed  the  uncertainties  were  less  than  5% 
in  laminar  region  and  10%  in  turbulent  region. 

On  the  basis  of  the  experimental  results  obtained  by  Bedecarrats 
et  al.  [123]  and  Grozdek  et  al.  [117],  Monteiro  and  Bansal  [124] 
evaluated  the  mentioned  four  rheological  models.  Among  these 
models,  the  Herschel-Bulkley  model  was  considered  as  the  best 
one  and  the  only  one  which  predicted  the  trend  of  experimental 
results  well  due  to  the  consideration  of  ice  particles  diameter. 


3.2.2.  MPCS 

It  is  believed  that  MPCS  can  be  treated  as  Newtonian  fluid  if 
the  slurry  is  homogeneous  due  to  the  very  fine  particles,  and  this 
conclusion  has  been  verified  by  many  measured  viscosities  using 
rheometer  [9-15,125-128],  as  given  in  Table  1,  except  Zhang 
et  al.  [129]  who  obtained  pseudo-plastic  feature  of  10-40  wt% 
MPCS  using  a  rotating  cylinder  type  of  rheometer. 

The  results  reveal  the  viscosity  of  MPCS  is  related  with  the 
particle  fraction  and  temperature.  Vand’s  equation  was  used  by 


many  researchers  [9,10,12,125,127]  to  predicted  the  viscosity  of 
MPCS,  given  by  the  following  equation: 

p  =  Pf(\-(p-Atp2)-25  (6) 

where  A  is  the  constant  that  concerning  the  particle  size,  shape 
and  rigidity,  which  is  1.16  with  the  assumption  that  the  particles 
are  ideal  rigid  sphere.  Dammel  and  Stephan  [125]  found  this 
equation  underestimated  the  viscosity  of  MPCS,  since  the  parti¬ 
cles  in  MPCS  was  obviously  not  rigid  sphere;  Yamagishi  et  al.  [9] 
used  A= 3.7  while  Wang  and  Lin  [127]  determined  the  value 
of  4.4  according  to  their  experimental  result.  Considering  as 
Newtonian  fluid,  the  calculated  friction  factor  using  Poiseuille 
equation  and  Blasius  equation  given  by  Yamagishi  et  al.  [9] 
showed  perfect  agreement  with  the  measured  result,  as  shown 
in  Fig.  17(a);  the  friction  factor  of  MPCS  in  turbulent  flow  region 
measured  by  Wang  et  al.  [12]  was  higher  than  the  predicted  value 
by  Blasius  equation,  as  shown  in  Fig.  17(b),  and  the  authors 
attributed  this  higher  value  to  the  high  tube  wall  roughness; 
Inaba  et  al.  [10]  reported  the  overestimation  of  Blasius  equation 
in  MPCS  turbulent  flow,  as  shown  in  Fig.  17(c),  the  authors 
explained  it  by  the  flow  drag  reduction  effect  of  the  fine  particles 
(1.5  pm)  and  found  that  the  increase  of  the  ratio  of  larger  particles 
(17  pm)  caused  the  smaller  discrepancy  between  the  measured 
friction  and  the  calculated  values  by  Blasius  equation. 

3.2.3.  CHS 

So  far,  CHS  was  determined  as  non-Newtonian  fluid  in  all  the 
open  literature  [24-33,36,37,131].  Based  on  the  measured  flow 
curves  of  TBAB  CHS  in  16  mm  tubes,  as  shown  in  Fig.  18(a), 
Darbouret  et  al.  [25]  believed  that  TBAB  CHS  behaved  as  Bingham 
fluid.  The  authors  presented  the  corresponding  yield  stress  and 
viscosity  which  were  about  0-14.18  Pa  and  0.008-0.17  Pa  s, 
respectively,  depending  on  the  volume  fraction  of  CHS,  and  found 
good  agreement  between  the  measured  pressure  drop  in  laminar 
flow  region  and  the  predicted  values  using  Buckingham-Reiner 
equation.  In  addition,  Song  et  al.  [27]  also  treated  TBAB  CHS  as 
Bingham  fluid  in  their  study  of  heat  transfer  characteristics  of  the 
slurry. 

Other  researchers  determined  TBAB  CHS  as  pseudo-plastic 
fluid  by  the  flow  test  in  pipes,  including  Hayashi  et  al.  [24],  Xiao 
et  al.  [26],  Ma  et  al.  [28],  and  Kumano  et  al.  [30].  The  results  of  the 
relation  between  ln(DAP/4L)  and  ln(8 v/D)  of  type  B  TBAB  CHS 
measured  by  Ma  et  al.  [28]  are  shown  in  Fig.  18(b),  which 
indicated  that  the  Power-law  model  was  appropriate  to  TBAB 
CHS.  The  flow  behavior  index  n  and  fluid  consistency  coefficient  K' 
obtained  by  different  authors  are  summarized  and  presented  in 
Fig.  19,  and  it  was  commonly  accepted  that  larger  slurry  fraction 
caused  the  slurry  deviated  Newtonian  fluid  more.  Using  the 
modified  Reynolds  number  (Eq.  (4)),  Ma  et  al.  [28]  calculated 
the  slurry  pressure  drop  in  straight  tubes:  Poiseuille  equation  was 
used  for  laminar  flow,  and  the  relative  error  was  +10%  and  +15% 
in  6  mm  and  14  mm  tube,  respectively;  while  Dodge-Metzner 
correlation  was  used  for  turbulent  flow,  and  the  relative  error  was 
-20%  to  -25%  and  +8%,  respectively. 

Based  on  the  tested  viscosities  using  rheometer,  Suzuki  et  al. 
[36,37]  obtained  the  shearing-thinning  pseudo-plastic  behavior  of 
TME  CHS,  and  Hashimoto  et  al.  [131]  obtained  the  similar 
behavior  of  TBAB  CHS  and  TBAF  CHS.  The  latter  study  also  found 
that  the  TBAB  and  TBAF  solution  were  both  Newtonian  fluid  when 
the  shear  rate  was  low  but  their  apparent  viscosity  increased 
when  shear  rate  was  larger  than  500  s  V 

3.3.  Pressure  drops  and  flow  friction  of  PCS 

Compared  to  single-phase  fluid,  pressure  drops  of  the  two- 
phase  PCS  are  dominated  by  many  other  parameters  in  addition 
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Table  1 

Measured  viscosities  of  MPCS. 


Authors 

Materials  (including 

PCM  material  and 
coating  material) 

Slurries 

Temperature, 

°C 

Viscosity,  mPa  s 

Descriptions 

Yamagishi 
et  al.  [130] 

n-Tetradecane; 
melamine  resin 
n-Tetradecane; 
melamine  resin;  1  wt% 
surfactant 

dp= 5  pm; 

<ps=20-30  vol% 

25 

3-1000  with  0.2-100  s-1  shear 
rate 

1-7 

p  decreased  with  increasing  shear  rate 

p  was  independent  of  shear  rate 

Yamagishi 
et  al.  [9] 

Octadecane;  melamine- 
formaldehyde 

dp— 2-10  pm; 

(ps= 7-30  vol% 

25 

1.17-10.4 

p  increased  with  particle  fraction;  nearly 
independent  of  temperature 

Alvarado 
et  al.  [11] 

n-Tetradecane;  gelatin 

dp=2-10  pm; 
xs= 4.6- 17.7  wt% 

2-12 

The  ratio  to  that  of  water  about 
1-3.5 

p  increased  with  particle  fraction,  and  was 
independent  of  temperature 

Wang 
et  al. 

[12,13] 

1  -Bromohexadecane 
amino  plastics 

dp  a=10.112  pm; 
xs= 5-27.6  wt% 

20 

1.57-8.42 

p  increased  with  particle  fraction 

Inaba  et  al. 
[101 

Small  particle; 
n-tetradecane 

Large  particle: 
n-docosane 

dp,a= 1.5  |xm; 
xs=20  wt%; 
dp,a=  17  pm; 
xs=0-50  wt% 
(compared  to 
small  particles) 

4-67 

0.5-5 

p  increased  with  particle  fraction  and  decreased 
with  temperature 

Chen  et  al. 
[14J 

Urea-formaldehyde 

dp =1-20  pm; 
xs=5-15.8  wt% 

10,  20 

1. 5-4.5 

Dammel 

and 

Stephan 

[125[ 

n- Eicosane; 

polymethylmethacrylate 

dp= 1 .5-12  pm; 
xs= 10-20  wt% 

30-60 

0.5-4 

Ho  et  al. 
[126[ 

n-Eicosane;  urea- 
formaldehyde 

dp=4-10  pm; 
xs=2-10  wt% 

10-50 

0.5-2. 1 

Wang  and 
Lin  [127] 

n- Octadecane; 

melamine- 

formaldehyde 

dp =0.3-3  pm; 
xs=5-20  wt% 

20-40 

1.5-6 

Zeng  et  al. 
[128[ 

1  -Bromohexadecane 
(C16H33Br);  amino 
plastics 

dp= 1-20  pm 
(8.2  pm  average); 
xs=5.0  wt%, 

10.0  wt%,  15.8  wt% 

14.3 

10  °C:  5.0  wt%  slurry: 

2.05,10.0  wt%  slurry:  2.23, 

15.8  wt%  slurry:  4.18 

p  decreased  with  the  temperature  rising 

Zhang 
et  al.  [129] 

Nonadecane  melamine 
resin 

dp =0-4.5  pm; 
xs= 10-40  wt% 

10-50 

10-30  wt%  slurry:  2-12;  40wt% 
slurry:  12-35 

p  increased  with  particle  fraction  and  decreased 
with  temperature;  however,  the  shear  stress  did  not 
vary  linearly  with  shear  rate,  therefore  the  fluid  was 
considered  as  pseudo-plastic  fluid 

Tumuluri 
et  al.  [15] 

Octadecane  gelatin 

MPCS  with  multi-walled 
carbon  nanotubes 

dp=2-10  pm;  xs=7, 
11  wt% 

11  wt%  MPCM; 

1.1  wt%  MWCNT 

17-32 

25-35 

7  wt%  slurry:  1. 8-2.2  times  as 
large  as  that  of  water;  1 1  wt% 
slurry:  2. 5-3.0  times  as  large  as 
that  of  water 

10-1000  with  0.2-100  s-1  shear 
rate 

p  was  almost  constant  with  respect  to  temperature; 
p  increased  with  particle  fraction;  MPCS  was 
Newtonian  fluid  up  to  18  wt%  fraction 

p  increased  with  temperature;  fluid  exhibited  non- 
Newtonian  behavior,  and  Herschel-Buckley  model 
was  used 

to  the  flow  velocity  and  tube  diameter,  including  particle  fraction, 
particle  diameter  and  shape  and  so  on.  Therefore,  the  former 
measured  results  of  PCS  pressure  drops  exhibit  irregularities.  The 
recent  experimental  results  of  PCS  pressure  drops  are  summar¬ 
ized  in  Tables  2-4. 


3.3.1.  Laminarizing  effect 

Most  of  the  experimental  results  showed  that  the  pressure 
drops  of  ice  slurry  were  higher  than  that  of  the  carrier  fluid  at  the 


same  flow  velocity  in  tubes  [108],  However,  Knodel  et  al.  [132] 
found  that  the  flow  friction  of  ice  slurry  with  2-3  mm  ice  crystals 
in  a  24  mm  stainless  steel  tube  reduced  with  increasing  ice  fraction 
in  a  nonlinear  manner.  The  friction  factor  was  close  to  the  value  of 
water  when  the  ice  fraction  was  low,  and  a  sudden  decline 
occurred  as  the  ice  fraction  increasing  from  2  to  4  wt%,  thereafter 
the  friction  factor  changed  to  another  constant  value.  The  authors 
characterized  the  decline  as  re-laminarization  phenomenon  by 
considering  the  flow  turned  from  turbulent  state  to  laminar  state 
as  ice  fraction  increasing.  This  phenomenon  can  be  explained  by 
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Fig.  17.  Friction  factor  prediction  of  MPCS,  (a)  Yamagishi  et  al.  [9],  10.1  mm  tube;  (b)  Wang  et  al.  [12],  4  mm  tube  and  (c)  Inaba  et  al.  [10],  15  mm  tube. 


the  laminarizing  effect  of  the  ice  crystals:  the  interactions  of  ice 
crystals  and  the  carrier  fluid  reduce  the  fluid  turbulence  level.  By 
comparing  the  results  with  their  former  experiments  of  ice  slurry 
in  larger  tube  diameter  (48  and  100  mm  tube  diameter,  8  wt% 
transition  ice  fraction),  they  believed  that  the  interactions  between 
ice  crystals  and  tube  wall  were  enhanced  when  the  tube  diameter 
was  small,  and  the  transitional  ice  fraction  was  low. 

The  re-laminarization  phenomenon  was  also  observed  when 
ice  slurry  flowing  through  10-20  mm  copper  tubes  by  Niezgoda- 
Zelasko  and  Zalewski  [119]  and  in  22  mm  inner  tube  of  double¬ 
tube  heat  exchanger  with  smooth  and  corrugated  surfaces  by 
Bedecarrats  et  al.  [123].  As  shown  in  Fig.  20,  ice  slurry  flowing  at  a 
certain  velocity  is  gradually  damped  from  fully  turbulent  to 
laminar  state  as  the  ice  fraction  increasing.  In  the  fully  turbulent 
region,  the  pressure  drop  is  nearly  constant  when  the  ice  fraction 
is  low,  and  then  gradually  increases  with  ice  fraction.  Bedecarrats 
et  al.  [123]  believed  that  the  liquid  phase  dominated  the  flow 
pressure  drop  when  the  ice  fraction  was  low,  therefore  the 
pressure  drop  did  not  change  significantly  in  the  turbulent  region. 
The  pressure  drop  increases  more  rapidly  in  the  laminar  region 
than  that  in  turbulent  region,  which  was  attributed  by  Bedecar¬ 
rats  et  al.  [123]  to  the  increasing  friction  of  ice  particles  on  the 
tube  wall.  Moreover,  pressure  drop  reduces  with  ice  fraction 
when  ice  slurry  flows  in  the  transitional  region,  and  the  transi¬ 
tional  ice  fraction  moves  towards  the  high  ice  fraction  with 
increasing  flow  velocity  as  clearly  shown  in  the  Fig.  20(a). 

Similar  phenomena  were  observed  in  the  flows  of  MPCS  and 
TBAB  CHS.  Yamagishi  et  al.  [9]  noticed  the  lower  pressure  drop  of 


30  vol%  MPCS  at  the  flow  velocity  of  2. 0-2. 5  m  s_1  compared  to 
that  of  water,  and  attributed  the  laminarizing  effect  to  the 
depression  of  longitudinal  eddy  diffusivity  of  slurry  since  the 
particle  size  was  smaller  than  the  turbulence  length.  The  authors 
also  pointed  the  laminarization  was  beneficial  for  pumping  power 
reduction.  Xiao  et  al.  [142]  and  Ma  et  al.  [28]  reported  the 
laminarization  during  the  flowing  of  TBAB  CHS  in  tubes,  the 
variations  of  pressure  drop  with  particle  fraction  at  the  certain  flow 
velocity  were  quite  similar  to  that  of  ice  slurry,  as  shown  in 
Fig.  21(a)  and  (b),  the  reduction  of  pressure  drop  were  both  noticed 
occurred  at  the  slurry  fraction  range  of  12-16  vol%.  Ma  et  al.  [28] 
also  noticed  that  the  laminarizing  effect  was  particle  fraction  and 
flow  velocity  dependent,  which  indicated  the  laminarizing  effect 
would  disappear  when  the  particle  fraction  and  flow  velocity  were 
both  high.  This  was  caused  by  the  enlarged  particle  interactions  and 
therefore  the  mechanical  friction  was  larger. 


3.3.2.  Influence  of  flow  passage  geometry  on  flow  friction 

Different  from  that  of  the  single  phase  fluid,  the  flow  friction  of 
solid-liquid  two-phase  slurry  is  composed  of  viscous  friction 
caused  by  the  fluid  viscosity  and  the  mechanical  friction  which 
is  resulted  from  the  interactions  between  solid  particle  and  the 
flow  boundary  [144].  Therefore,  the  flow  friction  of  PCS  is  largely 
related  to  both  the  slurry  fraction  and  flow  passage  geometry, 
such  as  passage  diameter,  roughness  and  pattern. 

Relationships  between  the  flow  friction  and  Reynolds  number 
of  ice  slurry  in  tubes  with  different  diameters  were  obtained  by 
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Fig.  18.  Rheological  feature  determination  of  TBAB  CHS,  (a)  wall  shear  stress  and 
shear  rate  of  TBAB  CHS  flowing  through  16  mm  tube  [25]  and  (b)  flow  curve  of 
TBAB  CHS  flowing  through  6  mm  tube  [28], 


Gozdek  et  al.  [117],  Niezgoda-Zelasko  and  Zalewski  [119]  and 
Ilian  and  Viedma  [133].  As  shown  in  Fig.  22(a)-(c),  the  experi¬ 
mental  results  obtained  by  these  studies  show  similar  friction 
factor-Reynolds  numbers  relationships  though  the  used  Reynolds 
numbers  were  defined  in  different  ways:  Grozdek  et  al.  [117] 
determined  the  Reynolds  number  using  the  viscosity  calculated 
by  Thomas  relation  (Refs),  Niezgoda-Zelasko  and  Zalewski  [119] 
used  the  Reynolds  number  based  on  the  Bingham  model  (ReB), 
while  Ilian  and  Viedma  [133]  employed  the  referenced  Reynolds 
number  ( Rec f)  corresponding  to  the  carrier  fluid  at  0  C.  In  the 
laminar  flow  region,  it  was  noticed  that  both  the  higher  ice 
fraction  and  the  larger  tube  diameter  resulted  in  larger  flow 
friction,  and  the  corresponding  deviation  from  that  of  the 
Newtonian  fluid  also  became  larger.  Grozdek  et  al.  [117]  illumi¬ 
nated  the  influence  of  the  tube  diameter  on  the  ice  slurry  flow:  ice 
slurry  flowing  through  smaller  tubes  tended  to  be  more  homo¬ 
geneous  at  larger  mean  velocity  but  still  in  laminar  state,  while 
the  mean  velocity  was  relatively  low  in  larger  tube  which  resulted 
in  heterogeneous  flow  with  more  ice  particles  at  the  tube  top  but 
less  at  the  bottom,  therefore  the  distorted  flow  pattern  led  to 
higher  flow  friction.  Kumano  et  al.  [120]  observed  the  similar 
influence  of  tube  diameter  on  the  flow  of  ice  slurry  while  narrow 
tubes  with  the  tube  diameter  of  4.3  mm,  7.5  mm  and  10.2  mm. 

Ma  et  al.  [28]  regarded  that  the  particle-tube  wall  interac¬ 
tions  was  enlarged  when  TBAB  CHS  flowing  through  tube  with 
smaller  diameter,  hence  the  fluidity  of  the  slurry  deviated  from 
Newtonian  fluid  much  more  compared  to  that  of  the  slurry 
flowing  in  larger  tube,  and  then  the  flow  friction  was  higher 
due  to  the  larger  mechanical  friction  loss  in  the  slurry.  In  their 
study,  the  measured  flow  friction  factor  in  6  mm  tube  was  about 
20-25%  larger  than  the  predicted  value  by  Dodge-Metzner 
correlation  as  mentioned  before,  while  the  prediction  error  in 
the  case  of  14  mm  tube  was  only  +8%  larger.  They  believed  this 
tube  wall  constraint  effect  on  slurry  flow  would  disappear  if  the 
tube  diameter  was  large  enough  compared  to  the  particle  size. 

Both  Bedecarrats  et  al.  [123]  and  Ilian  and  Viedma  [133] 
observed  reasonably  higher  pressure  drop  of  ice  slurry  in  corru¬ 
gated  tube  compared  to  that  in  smooth  tube.  In  the  former  study, 
different  from  the  slight  increase  in  smooth  tube,  the  pressure 


TBAB  CHS  fraction 


Ma  et  al.,  tpye  B  TBAB  CHS 
Copper  tube 
— ■  —  6  mm  tube 
— 14  mm  tube 


Xiao  et  al.,  tpye  B  TBAB  CHS 
21  mm  plastic  tube 
— 20  wt%  initial  solution 
— O—  33  wt%  initial  solution 


Kumano  et  al. 
Stainless  steel  tube 
tpye  A  TBAB  CHS 
—A—  4.3  mm  tube 
—A—  7.5  mm  tube 
—  10.2  mm  tube 


tpye  B  TBAB  CHS 
♦  4.3  mm  tube 
— O —  7.5  mm  tube 
<3>  10.2  mm  tube 


Fig.  19.  Summary  of  the  parameters  of  Power-law  model  for  TBAB  CHS,  including  researches  by  Xiao  et  al.  [26],  Ma  et  al.  [28]  and  Kumano  et  al.  [31]  (the  used  TBAB  CHS 
fraction  in  [26]  is  volume  fraction  while  in  [28,31]  are  mass  fraction). 


Table  2 

Pressure  drops  of  ice  slurries. 


Authors 

Base  liquids 

Slurries 

Geometries 

Flow  rates 

Correlations 

Descriptions 

Knodel  et  al. 
[132] 

Water 

dp =2-3  mm; 
xs  <  1 2  wt% 

Stainless  steel  tube, 

D,=24  mm,  L= 4.596  m 

2. 8-5.0  m  s~’ 

/  =  0.946  x  0.184Re-°2,  38,000  <  Re  <  74,000 

A P/L  was  nearly  constant  when  xs  <  2  wt%  and  xs  >  5  wt%; 

A P/L  reduced  when  2  <xs  <  5  wt%  due  to  the  laminarizing 
effect 

Meewisse 
and  infante 
Ferreira 
[134] 

5-39%  Various 
water-brine 

dp=0.1-l  mm; 
xs  <  30  wt% 

Circular  tube,  D,= 20  mm 

3  m  s_1 

A  P/L  increased  with  ice  fraction;  A  P/L  depended  on  the 
added  freezing  point  depressant.  Slurry  with  methanol, 
NaCl  and  potassium  formate  had  lowest  pressure  drop 

Bellas  et  al. 
[135] 

5%  Propylene 
solution 

dp = 0.1 7- 
0.6  mm; 
xs  <  20  wt% 

Plate  heat  exchanger,  24 
plates,  Dh= 4  mm, 

W=112  mm,  H=310  mm 

1.0-3. 7  m3  h~' 

1 

AP  increased  with  vs  and  xs;  only  a  3-8%  increment  in  AP 
was  observed  as  xs  varying  from  5.0  to  20.0  wt% 

Norgaard 
et  al.  [136] 

16%  Propylene 
glycol  solution 

xs  <  30  wt% 

Plate  heat  exchanger,  10  plates, 
Dh= 4  mm,  W=70  mm, 

H=300  mm 

50-300  Lh-' 

AP  of  ice  slurry  in  heat  exchange  generally  increased  with  vs 
and  xs;  slight  reduction  was  observed  when  flow  rate 
increased  from  100  to  150  L  h-1  with  inlet  30  wt%  ice 
fraction 

Hagg  [137] 

33.2  wt% 

HYCOOL 
solution  and 

35  wt%  ethanol 
solution 

xs  <  30  wt% 

Stainless  steel  tube,  D,= 9,  15, 
25  mm,  L=  1.2  m 

0.5-2  ms-1 

A  P/L  increased  with  vs  and  xs  and  with  reduced  pipe 
dimensions 

Niezgoda- 
Zelasko  and 
Zalewski 
[119] 

10.6%  Ethyl 

alcohol 

solution 

dp =0.05- 
0.35  mm, 
xs  <  30  wt% 

Copper  tubes,  D,=10,  16, 

20  mm,  L= 4.6  m 

0.1-4.5  m  s-1 

Buckingham-Reiner  equation  for  laminar  flow,  excellent 
agreement;  Tomita  method  for  turbulent  flow,  uncertainty 
was  better  than  10% 

A  P/L  increased  at  the  fully  turbulent  and  more  rapidly  at  the 
laminar  flow  regions  with  xs  at  a  certain  vs,  but  decreased  at 
the  transitional  region,  while  the  xtra  increased  with  vs:  Retra 
increased  with  xs 

Shire  et  al. 
[138] 

4.76  wt%  NaCl 
brine 

0.02-0.05  mm  in 
width,  0.1- 
0.2  mm  in  length, 
41  <  xs  <  67  wt% 

Borosilicate  glassware, 

D,= 25  mm,  L=1  m 

0.5-2  ms-1 

A  P/L  increased  with  vs  and  xs;  A  P/L  depended  largely  on  a 
critical  yield  stress  for  high  xs 

Ilian  and 

Viedma 

[139] 

9%  NaCl  brine 

dp  «  0.5  mm, 
xs  <  25  wt% 

Stainless  steel  corrugated  tube, 
D,=  16  mm,  L=1.25  m 

700, 

1200  kg  h_1 

/  =  6.12/c046Re^016,  k  is  the  severity  number  of  the 
corrugated  tube 

AP  was  almost  constant  with  xs 

Shire  et  al. 
[140] 

5%  NaCl  brine 

0.02-0.05  mm  in 
width,  0.1- 
0.2  mm  in  length, 

1 0  <  xs  <  60  vol% 

Plate  heat  exchanger,  12 
plates,  Dh= 7  mm,  W=250  m, 
H=1000  m,  tubular  heat 
exchanger,  D/,= 24  mm,  L=3  m 

0-2.5  Ls-' 

AP  increased  with  vs  and  xs;  AP  had  sharp  increase  when  xs 
was  higher  than  40  vol%  in  the  tubular  heat  exchanger; 
required  pressure  drop  was  lower  than  the  predicted  from 
models;  mixing  mitigated  the  block  caused  by  high  ice 
fraction  as  well  as  large  ice  particles 

Bedecarrats 
et  al.  [123] 

10%  Ethanol 
solution 

dp  «  0.1  mm, 
xs  <  30  wt% 

Smooth  and  corrugated  tube, 
Dj= 22  mm,  1=2  m 

OJ-LOms-1 

At  a  certain  vs,  A  P/L  increased  slowly  with  xs  within  smooth 
tube,  but  the  increment  stopped  at  critical  xs,  while  A  P/L 
was  nearly  constant  within  corrugated  tube;  then  A  P/L 
reduced  to  a  minimum  value  and  then  rapid  increase;  x^ 
increased  with  vs 

Grozdek 
et  al.  [117] 

10.3  wt% 

Alcohol 

solution 

dp  x  0.3  mm, 
xs  <  30  wt% 

Stainless  steel  tube,  D,= 9,  15, 
25  mm,  L=1.2  m 

0-2.6  ms-1 

The  authors  tried  several  kinds  of  correlations,  the  best  ones 
were  the  modified  Buckingham-Reiner  for  laminar  flow 
and  Dodge-Metzner  and  Tomita  for  turbulent  flow 

Ice  slurry  flow  behaviour  was  time  dependant;  A  P/L  generally 
increased  with  xs  at  a  certain  vs,  however,  the  delayed 
laminar-turbulent  transition  with  higher  xs  resulted  in  lower 
A  P/L  than  that  of  carrier  fluid  sometimes;  when  xs<  10- 
1 5  wt%,  the  slurry  could  be  treated  as  Newtonian  fluid 

9%  NaCl  brine 

dp  «  0.5  mm, 
xs  <  26  wt% 

Stainless  steel  tube, 
dJDi= 0.01 4-0.031 

0.2-5  ms-1 

Smooth  tube  (subscript  cf  stands  for  the  carrier  liquid  at 

0  °C):  Laminar  flow:  f  =  ARe°p 

/  increased  with  xs  in  laminar  region,  but  was  nearly 
constant  in  turbulent  flow  region;  /  increased  with  the 
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Table  2  ( continued ) 


Authors 

Base  liquids  Slurries  Geometries  Flow  rates 

Correlations 

Descriptions 

Ilian  and 

Viedma 

[133] 

A  =  0.00024(  1  -0.1xs/l  00)8’  8122 

T  /  .  x  0.6781  a 

x  3.95  exp((a  +  l)b) 

a  =  (0.1427-0.174  In  B), 
b  =  (0.22435-0.000308pc/)xs/100 

B  =  1.295D2pqfebtO°05’“p-/2  (Sr)”66. 
Turbulent  flow:  f  =  0.31 6cRe^0  25d, 
c  =  1 .4074(  1  -0.1xs/100)e-007185x- 
d  —  l.l-0.02847xs 

decrease  of  dm/D,  in  laminar  region,  but  was  nearly  constant 
in  turbulent  flow  region;  sharp  increase  was  found  when  xs 
was  up  to  26  wt%  in  corrugated  tube 

Kumano 
et  al.  [120] 

5  wt%  Ethanol  dp=0.1-0.4  mm,  D,= 4.3, 
solution  xs  <  25  wt% 

7.5,  10.2  mm,  L=2.5  m 

Homogeneous  flow  without  ice  particles  accumulated  at  the 
tube  top;  A  P/L  increased  with  xs  in  the  laminar  flow  region; 

A  P/L  was  constant  or  decreased  slightly  with  xs  in  the  turbulent 
flow  region  when  xs  was  low,  but  increased  with  xs  when  xs 
was  high;  slurry  flowing  through  larger  tube  deviated  more 
from  the  Newtonian  fluid  and  had  higher  flow  friction 

Fernandez- 
Seara  et  al. 
[118] 

10wt%  dp= 0.025-  Fan  coil  heat  exchanger,  0.3- 

Ethylene  glycol  0.25  mm,  smooth  copper  tube,  D,=3/8" 

solution  xs  <  20  wt% 

1.1  ms-1 

Poiseuille  model  for  laminar  flow  with  relative  errors  of  -  A P/L  increased  with  vs  and  xs 

21.5%  to  29%,  Petukhov  equation  for  turbulent  flow  with 
relative  errors  of  -6.3%  to  9%,  Blasius  equation  for  turbulent 
flow  with  relative  errors  of  -17%  to  1% 

Table  3 

Pressure  drops  of  MPCS. 

Authors 

PCM  and  shell  material 

Slurries 

Geometries 

Flow  rates 

Correlations 

Descriptions 

Yamagishi 
et  al.  [9] 

Octadecane;  melamine- 
formaldehyde 

dp =2-10  pm, 

(ps= 7-30  vol% 

Stainless  steel  pipe, 

D,=  10.1  mm,  L= 800  m 

0-3.0  ms-1 

Poiseuille  equation  for  laminar 
flow  and  Blasius  equation  for 
turbulent  flow  with  perfect 
agreement 

A  P/L  increased  with  vs;  A  P/L  of  7-20  vol%  slurry  was  higher  than  that  of  water 
and  increased  with  cps ;  laminarization  was  observed  when  <ps=30  vol%,  which 
caused  A  P/L  smaller  than  that  of  water  at  vs=2.0-2.5  m  s-1 

Rao  et  al. 
[141] 

n-Octadecane; 

polymethylmethacrylate 

dp= 1-5  pm, 
xs=5-20  wt% 

Rectangular  copper 
minichannels, 

2  mm  x  4.2  mm  x  1 50  mm 

0.1-1.4  m  s 

-i 

AP  increased  with  vs  and  xs 

Alvarado 
et  al.  [11] 

n-Tetradecane;  gelatin 

dp  =  2-10  pm, 
xs= 5.9- 15.2  wt% 

Copper  pipe,  D,=10.9 
(smooth),  8  (smooth), 

8  mm  (enhanced), 

L=12.2  m 

0.5-1.5  m  s_1  in 

10.9  mm  tube; 

1.2-2.5  m  s-1  in 

8  mm  tube 

A  P/L  increased  with  vs\  In  10.9  mm  smooth  tube,  A  P/L  was  slightly  higher 
than  that  of  water  (calculated  value);  In  8  mm  tube,  A  P/L  of  5.9  wt%  slurry 
was  smaller  than  that  of  water  (calculated  value)  due  to  the  turbulence 
suppress  effect  by  small  particles;  A  P/L  in  smooth  and  enhanced  tubes  were 
similar 

Chen  et  al. 
[14] 

1  -Bromohexadecane; 
urea-formaldehyde 

dp— \ -20  pm, 
xs=15.8  wt%  at 

10  and  20  °C 

Stainless  steel  pipe, 

D,= 4  mm,  L=1.85  m 

0.2-1. 0  m  s 

AP  increased  with  vs  and  increased  slightly  with  temperature;  AP  was 
higher  than  that  of  water  at  the  same  vs,  but  the  difference  reduced  with 
increasing  flow  velocity 

Dammel 

and 

Stephan 

[125] 

n- Eicosane; 

polymethylmethacrylate 

dp  =1.5-12  pm; 
xs= 10-20  wt% 

Rectangular  copper 
minichannels, 

2  mm  x  4.2  mm  x  390  mm 

0.2-0.8  kg  min-' 

AP  increased  with  vs  and  xs 
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Table  4 

Pressure  drops  of  CHS. 


Authors 

Initial 

solution 

concentration 

Slurries 

Geometries 

Flow  rates 

Correlations 

Descriptions 

Darbouret 
et  al.  [25] 

20  wt%, 

35  wt% 

TBAB  CHS,  Type 
A:  xs=4- 
55  wt%;  Type  B: 
xs= 28-43  wt% 

D,  =  16  mm,  L=1  m 

0-300  Lh-1 

Buckingham-Reiner 
equation  for  laminar 
flow 

AP  increased  with  xs,  and  increased  linearly 
with  flow  rate;  Bingham  fluid  was 
determined  for  TBAB  CHS 

Xiao  et  al. 
[26] 

22  wt%, 

30  wt% 

TBAB  CHS: 
cps=0- 
16.24  vol% 

PVC  pipe,  Di= 21  mm,  L= 2  m 

0-5.25  ms-1 

Laminar  flow : 
Poiseuille  equation; 
Turbulent  flow: 
fF=0.3226ReMR427' 

AP  increased  with  vs:  AP  sometimes 
decreased  with  <ps;  power-law  feature  was 
determined  for  TBAB  CHS 

Suzuki 
et  al. 

[36,37] 

25  wt% 

TME  CHS: 
xs=0-18.8  wt%; 
Ethoquad  0/12 
as  surfactant 

Dj= 5  mm,  L= 0.5  m 

Re  (calculated 
based  on  water 
viscosity)  was  in 
the  range  of 
2000-20,000 

Drag  reducing  by  the  used  surfactants  was 
effective  for  TME  CHS  when  the  flow  rate  and 
the  slurry  fraction  were  both  high 

Xiao  et  al. 
[142] 

Type  A  TBAB 
CHS:  xs=0- 
27  wt% 

PVC  pipe,  D,  =  16  mm,  [=2.8  m 

0-3.5  m  s-1 

AP  increased  with  vs ;  as  the  increase  of  xs, 
the  flow  of  CHS  was  in  the  turbulent  region 
first  and  AP  increased  slightly;  when  xs  was 
in  the  range  of  10-15  wt%,  laminarization 
was  observed,  and  AP  dropped  largely;  then 

AP  increased  rapidly  in  the  laminar  flow 
region 

Ma  et  al. 
[28] 

10  wt% 

Type  B  TBAB 
CHS:  (ps= 0- 
20  vol% 

Copper  tube,  D,= 6,  14  mm, 
1=1.5  m 

0-2.5  m  s-1  in 

6  mm  pipe,  0- 
1.2  m  s-1  in 

14  mm  tube 

Laminar  flow: 
Poiseuille  equation; 
Turbulent  flow: 
Dodge-Metzner 
correlation 

In  laminar  region,  AP  increased  slightly  with 
(ps  in  6  mm  tube  but  changed  little  in  14  mm 
tube;  in  the  transitional  and  turbulent 
regions,  the  laminarizing  effect  caused  a 
decrease  of  AP  with  cps,  however,  the 
laminarizing  effect  was  volume  fraction  and 
flow  velocity  dependent 

Kumano 
et  al.  [30] 

22.5  wt% 

8.5  wt% 

Type  A  TBAB 
CHS:  xs=0- 
20  wt% 

Type  B  TBAB 
CHS:  xs=0- 
20  wt% 

Stainless  steel  pipe,  D,= 4.3, 

7.5,  10.2  mm,  L=1  m 

AP  increased  with  vs;  AP  increased  with  xs  in 
laminar  flow  region  but  slightly  decreased  in 
turbulent  flow  region;  In  turbulent  flow 
region,  the  flow  characteristics  changed  at  a 
certain  flow  velocity  while  xs  was  about 

1 5  wt%,  then  AP  dropped  a  lot 

Hashimoto 
et  al.  [131] 

20  wt%  TBAB 
solution; 

18  wt%  TBAF 

solution 

TBAB  CHS: 
xs=0-70  wt%; 
TBAF  CHS: 
xs=0-42  wt%; 
Ethoquad  0/12 
as  surfactant 

The  addition  of  surfactant  caused  the 
drastically  reduction  of  the  viscosities  of  both 
TBAB  CHS  and  TBAF  CHS 

Ma  et  al. 
[33] 

10  wt% 

Type  B  TBAB 
CHS:  cps= 0- 
17.5  voI% 

Plate  heat  exchanger,  1 0 
plates,  D/,—4  mm,  W=  72  mm, 
H=310  mm 

2.5-13  L  min-1 

/=  51.002Re„)P, 

60  <  Re^R  <  400 

AP  increased  with  vs:  AP  was  always  higher 
than  water;  AP  was  almost  constant  when 

5  <  cps  <  13.6  vol%,  but  increased  a  lot  when 
q>s=  17.5  vol% 

Ma  et  al. 
129] 

15  wt% 

Type  B  TBAB 
CHS:  xs=0- 
30  wt% 

Copper  pipe  fittings,  14  mm 

90°  elbow,  14-6  mm  sudden 
contraction,  6-14  mm  sudden 
expansion,  14  mm  tee 

1-8  m  s-1  in 

6  mm  pipe, 
0.4-2.4  ms-1  in 
14  mm  tube 

The  local  loss  friction  increased  slightly  with 
xs  when  xs  was  low  but  increased  rapidly  if  xs 
was  larger  than  20  wt% 

Suziki  et  al. 
[143] 

35  wt% 

ammonium 

alum 

Ammonium 
alum  hydrate 
slurry: 
xs=12  wt% 

Copper  tube,  D,= 13  mm, 

L= 2  m 

0.5-1 .5  ms"1 

AP  of  slurry  was  slightly  larger  than  that  of 
water  at  the  same  Re.  surfactant  was 
effective,  with  which  AP  was  even  lower 
than  that  of  water 

drop  of  ice  slurry  flowing  in  the  turbulent  region  in  corrugated 
tube  was  nearly  constant  or  only  had  slight  reduction  as  increas¬ 
ing  the  ice  fraction.  In  the  latter  study,  pressure  drop  in  the 
laminar  region  was  noticed  a  sharp  increase  when  the  ice  fraction 
was  up  to  26  wt%,  and  the  authors  also  successfully  introduced 
the  severity  parameter  relevant  to  the  corrugated  geometrical 
characteristics  to  predict  the  pressure  drop.  Alvarado  et  al.  [11] 
compared  the  pressure  drops  of  5.9  wt%  MPCS  in  regular  and 
enhanced  8  mm  tubes  (the  tube  roughness  was  considered  as 
0.1  mm),  however,  the  results  showed  the  similar  pressure  drop 
in  magnitude. 


3.4.  Heat  transfer  characteristics  of  PCS 

Numerous  experiments  have  been  performed  to  investigate 
the  overall  and  local  heat  transfer  coefficients  of  PCS  flowing 
through  the  heated  tube  and  heat  exchangers,  the  relating  investi¬ 
gations  were  summarized  in  Tables  5-7.  It  was  expected  that  the 
latent  heat  of  the  solid  particles  in  slurry  could  improve  the  heat 
transfer  of  ice  slurry,  and  many  reported  investigations  validated 
this  improvement,  as  shown  in  Fig.  23  which  presented  the  heat 
transfer  enhancement  of  ice  slurry  summarized  by  Stamatiou  and 
Kawaji  [112],  however,  the  comparison  may  conducted  under  the 
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Fig.  20.  Pressure  drops  of  ice  slurry  as  the  increase  of  ice  fraction,  (a)  16  mm 
copper  tubes  [119]  and  (b)  22  mm  inner  tube  of  double-tube  heat  exchanger  with 
smooth  surface  [123], 


same  Reynolds  numbers  since  the  flow  turbulence  of  slurry  was 
damped  by  the  particle  flow.  Some  researchers  also  attributed 
this  heat  transfer  coefficient  improvement  to  the  turbulence  in 
the  near-wall  thermal  boundary  layer  due  to  the  movement  of 
solid  particles.  However,  it  was  believed  that  particle  fraction, 
flow  velocity,  heat  flux,  particle  shape  and  size,  pipe  shape  and 
dimension,  heat  exchanger  type  and  size  imposed  influences  on 
the  heat  transfer  of  ice  slurry,  therefore  the  combined  effect  can 
be  complicated. 


3.4.1.  Local  heat  transfer  along  the  flow  direction 

Stamatiou  and  Kawaji  [112]  obtained  the  variations  of  local 
Nusselt  numbers  of  ice  slurry  flowing  though  heated  rectangular 
channels,  as  shown  in  Fig.  24(a).  Nuloca/  was  large  at  the  channel 
inlet  and  decreased  slightly  towards  the  exit  of  the  channel.  The 
authors  attributed  the  large  Nu(oca(  at  the  inlet  to  the  developing 
flow  and  the  following  decrease  to  the  more  developed  flow.  Hagg 
[137]  reported  the  similar  slight  reduction  of  local  heat  transfer 
coefficient  of  ice  slurry  along  the  flow  direction  in  21  mm 
stainless  steel  pipe,  as  shown  in  Fig.  24(b).  The  authors  believed 
that  the  phase  change  and  turbulence  caused  by  ice  particles  in 
the  thermal  boundary  layer  led  to  smaller  temperature  gradient 
at  the  entrance  and  therefore  high  hiocai,  and  the  turbulence 
reduced  as  the  ice  particles  melting  along  the  flow  direction  and 
therefore  hiocai  decreased.  Moreover,  lonescu  et  al.  [147]  reported 
hiocai  of  ice  slurry  decreased  greatly  at  the  channel  inlet,  followed 
by  nearly  constant  values. 


Fig.  21.  Pressure  drops  of  TBAB  CHS  as  the  increase  of  particle  fraction,  (a)  type  A 
TBAB  CHS  in  27  mm  PVC  tube  [142]  and  (b)  type  B  TBAB  CHS  in  6  mm  copper  tube 
[28]. 


Yamagishi  et  al.  [9],  Alvarado  et  al.  [11],  Wang  et  al.  [12,13] 
and  Tumuluri  et  al.  [15]  obtained  similar  hiocal  variations  of 
turbulent  flow  of  MPCS  in  heated  tubes,  typical  results  were 
shown  in  Fig.  25.  The  local  heat  transfer  coefficient  increased  at 
the  tube  inlet,  where  the  slurry  temperature  was  lower  than  its 
phase  change  temperature  (supercooled  slurry).  The  increasing 
hiocai  was  caused  by  the  gradually  increasing  melted  particles  near 
the  tube  wall;  then  the  slurry  temperature  reached  its  phase 
change  temperature,  and  hiocai  decreased  in  this  region  since  the 
number  of  solid  particles  decreased  gradually  and  the  fully 
melted  layer  near  the  tube  wall  deteriorated  the  heat  transfer 
from  tube  wall  to  fluid  core.  After  the  slurry  was  finally  fully 
melted,  hiocai  increased  again  due  to  the  rising  slurry  temperature. 
Yamagishi  et  al.  [9]  and  Wang  et  al.  [12,13]  also  showed  the 
variation  of  /i/oca(  in  laminar  flow  region,  which  generally 
decreased  along  the  flow  direction.  The  former  study  pointed 
that  the  phase  change  particles  in  the  laminar  flow  was  possibly 
not  completely  melted,  leading  to  generally  decreasing  solid 
particle  number  and  then  decreasing  heat  transfer  coefficient. 

As  shown  in  Fig.  26(a),  Kumano  et  al.  [31]  reported  the 
generally  decreasing  h/ocai  of  TBAB  CHS  in  laminar  flow  region, 
and  attributed  the  decrease  to  the  undeveloped  boundary  layer  of 
the  slurry.  The  authors  also  found  that  the  variation  of  htocal  along 
the  heating  tube  changed  little  under  the  turbulent  flow  condi¬ 
tion,  which  was  caused  by  the  rapidly  developed  boundary  layer. 
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Fig.  22.  Relationships  between  Reynolds  number  and  friction  factor  of  ice  slurries,  (al )  and  (a2)  Niezgoda-Zelasko  and  Zalewski  (119]  and  (b)  Grozdek  et  al.  [1 17J;  (cl)  and 
(c2)  Ilian  and  Viedma  [133], 


The  results  reported  by  Ma  et  al.  [28]  were  shown  in  Fig.  26(b), 
and  it  can  be  seen  from  the  figure  that  /iloca/  decreased  along  the 
heating  tube  but  increased  after  all  the  solid  particles  were  fully 
melted.  This  variation  was  explained  by  the  similar  reason  to  that 
of  MPCS,  except  there  was  no  supercooling  region  for  TBAB  CHS. 


3.4.2.  Influence  of  particle  fraction  on  heat  transfer 

Heat  transfer  performances  of  ice  slurry  flowing  through 
heated  circular  tubes  have  been  investigated  by  many  research¬ 
ers.  However,  sometimes  quite  different  even  opposite  character¬ 
istics  of  the  heat  transfer  coefficients  with  ice  fraction  were 
reported,  and  the  same  phenomena  happened  to  MPCS  and  CHS. 

Similar  to  the  re-laminarization  phenomenon  of  the  ice  slurry 
flow  as  aforementioned,  Knodel  et  al.  [132]  also  observed  the 
reduction  of  the  Nusselt  number  of  ice  slurry  in  the  heated  tube, 
as  shown  in  Fig.  27(a).  This  reduction  was  attributed  to  the 


damped  flow  turbulent  degree  by  laminarzing  effect.  As  the 
reported  by  Niezgoda-Zelasko  [145]  and  Niezgoda-Zelasko  and 
Zelasko  [148],  there  was  a  slight  increment  of  the  heat  transfer 
coefficients  of  ice  slurry  as  ice  fraction  increasing  from  5  to 
30  wt%  in  the  laminar  flow  region,  while  the  heat  transfer 
coefficients  were  nearly  constant  in  the  turbulent  flow  region, 
as  shown  in  Fig.  27(b).  Clearly,  there  is  a  sharply  decline  of  the 
heat  transfer  coefficients  as  ice  slurry  flowing  in  the  transitional 
region.  Moreover,  the  authors  pointed  that  the  influence  of  ice 
fraction  on  heat  transfer  was  more  distinct  in  the  laminar  flow 
region  while  the  heat  transfer  coefficients  were  more  strongly 
dominated  by  flow  velocity  than  by  ice  fraction  in  the  turbulent 
flow  region.  The  same  conclusion  was  also  obtained  by  Grozdek 
et  al.  [149]  and  Kumano  et  al.  [150]  in  the  case  of  heated  tubes, 
and  Bedecarrats  et  al.  [123]  in  the  case  of  double-tube  heat 
exchanger.  Grozdek  et  al.  [149]  observed  generally  increased  heat 
transfer  coefficient  with  ice  fraction  as  well  as  flow  rate,  as  shown 


Table  5 

Heat  transfer  of  ice  slurries. 


Authors  Base  liquids  Slurries  Geometries 


Knodel  et  al. 
[132] 

Water 

dp =2-3  mm,  Stainless  steel  tube,  D,= 24  mm, 

xs<12wt%  L=4.596  m 

Bellas  et  al. 
[135] 

5%  Propylene  solution 

dp =0.1 7-0.6  mm,  Plate  heat  exchange,  24  plates, 
xs  <  22  wt%  D/,= 4  mm,  W=  112  mm,  H=310  mm 

Norgaard  et  al. 
[136] 

16%  Propylene  glycol 
solution 

xs  <  30  wt%  Plate  heat  exchange,  10  plates, 

Dh= 4  mm,  W=70  mm,  H= 300  mm 

Stamatiou  and 
Kawaji  [112] 

6.2  wt%  NaCl  brine 

dp =0.1  -0.2  mm,  Rectangular  channels, 

xs  <  25  wt%  533  mm  x  305  mm  x  25  mm, 

533  mm  x  100  mm  x  13  mm 

Niezgoda- 
Zelasko  [145] 

10.6%  Ethanol  solution 

dp=0.1-0.15  mm,  Copper  tubes,  D,=10,  16,  20  mm, 
xs<30wt%  L=4.6  m 

Lee  et  al.  [146] 

6.5%  Ethylene  glycol 
solution 

dp  ss  0.27  mm, 
xs  <  26  wt% 

Double-tube  heat  exchanger,  inner 
tube:  D,= 13.84  mm,  outer  tube: 
D,=25  mm,  L=1.5  m 

Ionescu  et  al. 
[147] 

Composed  by  ice 
particles  and  oil, 
dp  x  0.1  mm, 
xs=20,  25,  30  wt% 

Rectangular  channel, 

8  mm  x  80  mm  x  1000  mm 

Flow  rates 
2.8-5.0  m  s-1 

1. 0-3.7  m3  h~ 

50-300  Lh-1 


0.1 -4.5  ms 

800- 

3500  kg  m-2 


Correlations 

Comments 

1  Num  =  0.885NuPetukhov 

Num  reduced  when  xs  <  6-8  wt% 

-1 

Laminar  flow : 

Nu(x)  =  4.0Czc/0486x°3(/itf  /nw)°  24 

hoveraii  remained  fairly  constant  as  the  increase  of  xs;  hoveraii 
increased  with  vs;  30%  increment  in  heat  transfer  capacity 
was  observed  as  the  xs  varying  from  5.0  to  20.0  wt% 

hoveraii  increased  with  vs  and  xs;  increment  with  xs  became 
little  at  high  vs 

Nuiocai  was  high  at  the  inlet  and  decreased  slowly  along  the 
flow  direction;  Nuiocai  and  Num  increased  with  xs  and  heat 
flux 

Turbulent  flow:  Num  =  NiiCn[l  +1.85x 
105x?72Re^-3(^//iw)2-47] 

Nucn  was  the  single-phase  Nusselt 
number  calculated  by  Gnielinski 
equation 

1  Laminar  flow:  Nu  =  0.78(Re  Pr  <5)° 72 

1  AH  Axy0027  (dpy 029 
\CpAT100y  VA/ 

h  changed  little  with  heat  flux  density;  h  increased  slightly 
with  xs  in  laminar  region;  h  decreased  with  xs  at  the 
transitional  region;  h  was  more  strongly  dependent  on  vs 
than  on  xs  in  turbulent  region;  h  of  10-20  wt%  ice  slurry 
can  be  lower  than  that  of  single-phase  carrier  fluid;  Nu  was 
much  larger  than  that  of  carrier  fluid 

Turbulent  flow:  Nu  =  0.031  (Re  Pr)0  73 
(  AH  Axy0018(dp\021 
\CP  AT  i  ooy  [d  J 

3  was  the  Bingham  term 

s-> 

At  low  flow  rate,  hoveraii  did  not  change  with  xs  when  xs  was 
smaller  than  about  10%,  then  increased  sharply  with  xs; 
hoveraii  increased  gradually  with  xs  at  high  flow  rate  range; 
hot  water  temperature  had  little  impact  on  hoverau 

Num(x )  indicated  Nu(x)  integrated 
from  0  to  x 

h local  had  a  great  decrease  at  the  inlet,  followed  by  a  nearly 
constant  value;  /i/oca/  increased  with  vs  and  xs 

Laminar  flow : 

-for  a  constant  wall  temperature: 
Num(x)  =  0.34  12Re042  Pr033 

(1  +xs)23 

-for  a  constant  heat  flux: 

Num(x>  =  0.22  149  Re0  51  Pr033 

<l+xs)27 

Transition  flow: 

Num(x)  =  0.03  35Re° 75  Pr° 33 

(l+xs)2-45 
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Niezgoda- 
Zelasko  and 
Zelasko  [148] 

10.6%  Ethanol  solution 

dp =0.1  -0.1 5  mm, 
xs  <  30  wt% 

Copper  tubes,  D,=10,  16,  20  mm, 
L=4.6  m,  rectangular  tube, 

0.0078  m  x  0.0265  m  x  3  m, 
rectangular  (slit)  tube, 

0.003  m  x  0.0358  m  x  2  m 

0.1-4.5mr' 

Gozdek  et  al. 
[149] 

10.3%  Alcohol  solution 

dp  k  0.3  mm, 
xs  <  22  wt% 

Stainless  steel  tube,  D,=21  mm, 

L= 0.94  m 

0-1.25  ms-1 

Bedecarrats 
et  al.  [123] 

10%  Ethanol  solution 

dp  «  0.1  mm, 
xs  <  30  wt% 

Double-tube  heat  exchanger  with 
smooth  and  corrugated  tubes,  inner 
tube:  Dj= 22  mm,  L=2  m,  outer 
tube:  D,= 33  mm,  L=1.5  m 

0.3-1. 9  ms"1 

Ilian  and 
Viedma  [158] 

9%  NaCl  brine 

dp  «  0.5  mm, 
xs  <  26  wt% 

Fan-coil  heat  exchanger,  cooper  tube, 
D,= 8.8  mm 

1000-4600  kg 

h-1 

Ilian  and 
Viedma  [133] 

9%  NaCl  brine 

dp  «  0.5  mm, 
xs  <  26  wt% 

Stainless  steel  tube,  dm/D, =0.014- 
0.031 

0.2-5  m  s~' 

Fernandez- 
Seara  et  al. 
[118] 

10  wt%  Ethylene  glycol 
solution 

dp =0.025- 
0.25  mm, 
xs  <  20  wt% 

Fan  coil  heat  exchanger,  smooth 
copper  tube,  D,=3/8" 

400,  600, 

800  kg  h-1 

Kumano  et  al. 
[150] 

5  wt%  Ethanol  solution 

dp  «  0.15  mm, 
xs  <  25  wt% 

D,=4.3,  7.5,  10.2  mm,  L= 2.5  m 

Turbulent  flow: 

-for  a  constant  wall  temperature: 
Num  (x)  =  (1 . 57  +  26.  lx°  75 ) 

£z0.33  +  1.76x2-1.1xs 

-for  a  constant  heat  flux: 

Num(x)  =  (2.06  +  78.8x]  85) 

£Z0.33  +  0.55x2  _o.54xs 


Laminar  flow: 
Nu=A&m(t$rm) 


WOV), 


y 


Similar  to  Niezgoda-Zelasko  [145];  smaller  hydraulic 
diameter  of  channel  caused  smaller  increase  in  h  of  ice 
slurry  compared  with  h  of  carrier  fluid 


Turbulent  flow: 

Nu=A(Re  Pr)m^y All  parameters 
are  seen  in  [148] 

Laminar  flow: 

Nu(x)  =  87.3  Gz° 1  (xs/l  00)°  67 


Turbulent  flow: 

^  =  1  +  0.001 03xs —2.003 

j-j30-xs/100 


n  1  rpju-xs/iuu  0.339 Recf 

Re"  30  (Xs/100)-^ 


h  changed  little  with  heat  flux;  h  generally 
increased  with  xs  and  vs ;  xs  had  greater  impact  on  h  than 
vs ;  h  of  10-15  wt%  ice  slurry  had  little  increase  or  no 
increase  in  comparison  with  carrier  fluid;  heterogeneous 
flow  pattern  of  ice  slurry  in  laminar  region  caused  higher  h 


At  a  certain  vs  and  as  increasing  xs,  hoveraii  increased  up  to  a 
critical  xs  within  smooth  tube,  followed  by  a  reduction  to  a 
minimum  value  and  then  a  rapid  increase;  hoveran 
decreased  first  within  corrugated  tube,  then  increased; 
hoveraii  depended  on  heat  flux  and  flow  pattern;  hoveran  was 
more  strongly  dependent  on  vs  than  xs  in  turbulent  region 

Heat  transfer  rate  increased  with  vs;  frost  formation, 
partial  blockage  and  excessive  diameter  ratio  caused  the 
worse  behaviour  prediction 

Smooth  tube:  Num  decreased  with  xs 

Nu  =  (0.845 -0.01 2 16xs  )Nuc/, 

Nucf  =  0.0344(Recf- 1500)°  78  Pr°cf7 
Corrugated  tube  ( k  was  the  corrugated 
parameter): 

Nu  =  (0.892-0.01 824xs  )Nucf, 

Nucf  =  0.374K°25(Rec/-1500)°  74Pr°cf4 

Heat  transfer  rate  of  ice  slurry  was  higher  than  chilled 
water  by  the  factor  of  3.7-4.6  depending  on  ice  fraction 
and  fan  rotation  rate;  80-88%  of  the  overall  thermal 
resistance  was  on  the  air  side 

Laminar  flow:  h  increased  with  xs  and  vs ;  h  depended  more  strongly  on  xs 

Nu  —  50  IRe00741  (  Xs  l0822 foA0292  than  on  Vs  *n  laminar  flow  region;  in  turbulent  region,  h 
~  M  ^100'  \^/  was  nearly  constant  at  low  xs  but  increased  largely  with  vs; 

h local  changed  little  with  distance 


Turbulent  flow:  Nu  =  0. 023 Re  m°  8  Pr'/3 
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Table  6 

Heat  transfer  of  MPCS. 


Authors 

Materials 

Slurries 

Geometries 

Flow  rates 

Correlations 

Comments 

Goel  et  al. 
[151] 

n- Eicosane,  urea- 
formaldehyde 

dp=  100,  250  pm, 

(ps= 5-20  vol% 

Copper  tube,  D,=3.14  mm, 

L=0.4m 

Reynolds  number 
was  200,  1000 
(only  for  5  vol% 
slurry) 

Stefan  number  was  the  most  dominant  factor  on  heat 
transfer;  lower  Stefan  number  (or  heat  flux)  caused 
lower  wall  temperature;  effect  of  cps  for  a  given  Stefan 
number  could  be  expected  when  cps  was  5  vol%;  larger 
particle  led  to  lower  wall  temperature;  the  effect  of 
homogeneity  could  be  negligible 

Yamagishi 
et  al.  [9] 

Octadecane;  melamine- 
formaldehyde 

dp=2-10  pm, 

(ps= 7-30  vol% 

Stainless  steel  tube,  D,=10.1  mm, 
L=800  D, 

0-3.1  ms-1 

In  laminar  flow  region,  h  decreased  along  the  flow 
direction;  in  turbulent  flow  region,  h  increased  in  the 
supercooling  region  along  the  flow  direction,  and  then 
decreased  in  the  melting  region,  then  increased  after 
finishing  the  melt;  higher  heating  power  led  to  lower  h\ 
h  generally  increased  with  xs  and  vs 

Inaba  et  al. 
[10] 

Small  particle:  n- 
Tetradecane;  large 
particle:  n-Docosane 

dp.a= 1-5  pm, 
xs=20  wt%; 
dp=  17  pm,  xs=0- 
50  wt %  (compared  to 
small  particles) 

Stainless  steel  tube,  D,=15  mm, 

L— 5.85  m 

300  <  Reyvw?  <  8000 

Nuiocai  decreased  along  the  flow  direction  and  reached 
nearly  constant  when  fully  melted;  Num  increased  with 
vs  and  the  ratio  of  large  particle 

Alvarado 
et  al.  [11] 

n-Tetradecane;  gelatin 

dp =2- 10  pm, 
xs= 5.9-1 5.2  wt% 

Copper  tube,  D,=10.9  (smooth),  8 
(smooth),  8  mm 
(enhanced)L=12.2  m 

The  variation  of  h  along  the  flow  direction  was  similar  to 
the  result  in  [6];  h  increased  significantly  with  vs ; 
increase  of  xs  caused  lower  flow  turbulence  and  then 
lower  h,  and  h  was  lower  than  that  of  water  at  the  same 
flow  velocity;  the  enhanced  tube  did  improve  the  heat 
transfer,  however,  the  enhancement  effect  was  weak  in 
the  case  of  high  xs 

Wang 
et  al. 

[12,13] 

1  -Bromohexadecane, 
amino  plastics 

dp>a= 10.112  pm, 
xs= 5-27.6  wt% 

Stainless  steel  tube,  D,= 4  mm, 

L= 1.46  m 

400  <  Rem  <  3500 

Laminar  flow:  Nu,ocai  =  cNusha„  a„d  London, 
c  was  1.336-1.418  in  the  case  of  5-15.8  wt% 
MPCS,  +15%  relative  error 

Num  =  0.8148 Reg,4593  Pr°m4836 

Ste- 01277[(Ll  +l2)/D]0.3059  was  the 

phase  change  region,  +10%  relative  error 
for  5-27.6  wt%  MPCS 

Turbulent  flow:  Num  =4.8527  x 

1 0  4Re{J;7733  Pr2m™'  Ste03159 

](ii  +L2)/Dr0333(pm/pwr24349 

+ 10%  relative  error  for  5-10  wt%  MPCS 

In  laminar  flow  region,  h  decreased  along  the  flow 
direction  and  reached  a  nearly  constant  value  in  the  fully 
melted  region;  Nu/oca/  increased  with  xs  and  was  almost 
independent  of  heating  power;  In  turbulent  flow  region, 
h  decreased  along  the  flow  direction  in  the  supercooling 
region,  and  then  decreased  in  the  melting  region,  then 
increased  in  the  fully  melted  region;  higher  heating 
power  caused  higher  h  in  the  last  region;  Num  of  5  wt% 
was  higher  than  that  of  10  wt%  MPCS  at  same  Re 

Rao  et  al. 
[152] 

n- Octadecane, 
polymethylmethacrylate 

dPta= 4.97  pm, 
xs=5-20  wt%; 

Rectangular  copper  minichannels, 

2  mm  x  4.2  mm  x  150  mm 

0.05- 

0.35  kg  min-1 

hm  increased  with  vs ;  hm  increased  with  xs  when  vs  was 
0.05  kg  min-1,  but  decreased  with  xs  when  vs  was  high 

Diaconu 
et  al. 
[153,154] 

Rubitherm  RT6, 
polycyclic  cell 

dp^ 2.24  pm 

Helically  coiled  copper  tube  in 
tank,  D,= 10  mm,  140  mm  helix 
diameter:  Water  flowed  inside  tube 
while  MPCS  was  contained  in  tank 

Nature  convection  correlation, 

Num  =  0.802Ra°278,  10n<Ra<1013 

Wang  and 
Lin  [127] 

n-Octadecane, 

melamine- 

formaldehyde 

dp =0.3-3  pm, 
xs=5-20  wt% 

Aluminum  rectangular  tube, 

12.8  mm  x  1.8  mm  x  900.2  mm 

4.17-10.0  gs-1 

Based  on  the  dimensionless  wall  temperature,  MPCS  had 
better  heat  transfer  performance  than  that  of  water 
except  at  very  small  dimensionless  axial  distance 

Wang 
et  al.  [155] 

dp =0.5-5  pm, 
xs=5-20  wt%,  with 
0.5  wt%  Ti02 
nanoparticles 

Copper  tube,  D,= 3.97  mm, 

1=1.834  m 

Nu local  decreased  along  the  flow  direction;  NuiOCai  was 
close  to  that  of  water;  the  heat  transfer  was  enhanced  by 
Ti02  nanoparticles,  and  the  enhancement  ratio  increased 
with  xs 
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in  Fig.  27(c).  There  were  ice  fraction  thresholds  at  about  10- 
15wt%,  after  which  the  heat  transfer  coefficient  increased  dras¬ 
tically.  The  authors  believed  that  the  transition  of  ice  slurry  flow 
pattern  from  homogeneous  to  heterogeneous  caused  the  accu¬ 
mulating  of  ice  particles  near  the  top  tube  wall  and  the  direct 
contact  of  ice  particles  with  tube  wall  enhanced  the  heat  transfer. 
On  the  other  side,  Ilian  and  Viedma  [133]  observed  that  the 
Nusselt  number  decreased  with  the  ice  fraction  at  the  same 
Reynolds  number,  as  shown  in  Fig.  27(d),  and  same  conclusions 
were  reached  when  using  corrugated  tube. 

The  results  given  by  Yamagishi  et  al.  [9]  showed  hiocai  of  MPCS 
flowing  through  heated  tube  was  higher  than  that  of  pure  water 
and  increased  with  the  slurry  fraction,  however,  the  comparison 
should  be  conducted  at  the  same  Reynolds  number.  The  flow 
turbulence  of  slurry  would  be  damped  as  increasing  particle 
fraction,  which  caused  the  reduction  of  heat  transfer  coefficient, 
therefore  /i/oca/  of  MPCS  decreased  as  increasing  particle  fraction 
under  the  condition  of  same  flow  velocity.  Ma  et  al.  [28]  noticed 
the  same  phenomena  when  studying  TBAB  CHS  flowing  through 
the  heated  tubes.  However,  other  researchers  reported  different 
results  about  the  heat  transfer  of  MPCS  in  the  heated  tubes:  Wang 
et  al.  [12,13]  found  the  Nusselt  number  of  5  wt%  MPCS  was  higher 
than  that  of  10wt%  MPCS  even  at  the  same  Reynolds  number 
(2100-3452);  Rao  et  al.  [152]  reported  the  mean  heat  transfer 
coefficients  of  MPCS  in  a  minichannel,  which  increased  as  slurry 
fraction  at  low  flow  velocity  (0.05  kg  min-1)  but  decreased  in  the 
case  of  high  flow  velocity  (0.15-0.35  kg  min-1);  Tumuluri  et  al. 
[15]  measured  the  heat  transfer  coefficient  of  7  wt%  and  11  wt% 
MPCS  in  heated  tube,  showing  /itoca/  of  1 1  wt%  MPCS  was  nearly 
the  same  as  that  of  7  wt%  MPCS  if  the  Reynolds  number  was  the 
same.  However,  fi(oca(  decreased  distinctly  with  slurry  fraction  at 
the  same  flow  velocity. 

Song  et  al.  [27]  obtained  quite  similar  variation  of  heat  transfer 
coefficient  of  type  A  TBAB  CHS  flowing  in  heated  tube  to  that  of 
ice  slurry  reported  by  Niezgoda-Zelasko  [145],  as  comparing 
Fig.  28(a)  to  27(b).  At  a  certain  flow  velocity  and  as  the  increase 
of  slurry  fraction,  the  flow  of  TBAB  CHS  was  gradually  damped 
from  turbulent  state  to  laminar  state,  therefore  the  heat  transfer 
coefficient  changed  little  with  low  slurry  fraction  (xs<10wt%) 
since  the  flow  velocity  dominated  the  heat  transfer  performance 
in  the  turbulent  flow  region.  The  re-laminarization  occurred  at 
10  wt%  <  xs  <  1 5  wt%  which  caused  large  reduction  of  heat  trans¬ 
fer  coefficient,  thereafter  the  heat  transfer  coefficient  increased  in 
the  laminar  flow  region.  Kumano  et  al.  [31]  also  reported  the 
similar  results  of  TBAB  CHS,  as  shown  in  Fig.  28(bl)  and  (b2), 
however,  the  re-laminarization  occurred  at  the  slurry  fraction  of 
about  15-17  wt%. 


3.4.3.  Influence  of  flow  passage  geometry  on  heat  transfer 

Niezgoda-Zelasko  and  Zelasko  [148]  discussed  the  impacts  of 
tube  diameter  on  the  heat  transfer  coefficients  based  on  the  test 
results  of  ice  slurry  flowing  through  circular  tubes  and  rectan¬ 
gular  channels.  Smaller  hydraulic  diameter  enhanced  the  micro¬ 
convection  effect  of  ice  particles  and  therefore  improved  the  heat 
transfer  coefficient.  However,  other  factors  bring  disadvantages 
for  the  heat  transfer  in  the  case  of  smaller  tube:  Reynolds  number 
at  the  same  flow  velocity  was  lower  than  that  in  the  case  of  larger 
tube,  and  the  melted  near-wall  layer  was  thicker  which  hindered 
the  heat  transfer  between  tube  wall  and  the  fluid  core.  Thus,  the 
authors  believed  that  the  smaller  hydraulic  diameter  the  smaller 
increase  of  heat  transfer  coefficients  compared  to  the  single¬ 
phase  carrier  fluid.  Results  reported  by  Kumano  et  al.  [150]  about 
ice  slurry  also  supported  this  conclusion. 

The  application  of  corrugated  tube  by  Bedecarrats  et  al. 
[123]  enhanced  the  heat  transfer  of  the  flowing  ice  slurry,  about 
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Table  7 

Heat  transfer  of  CHS. 


Authors 

Initial 

solution 

concentration 

Slurries 

Geometries 

Flow  rates 

Correlations 

Comments 

Takao 

TBAB  CHS 

Plate-fin  heat 

1-11  kg  min-1 

hm  generally  increased  with  vs ;  hm  was 

et  al. 

with  cold 

exchanger, 

larger  than  that  of  water  at  1-6  kg  min-1 

[21] 

carry  capacity 

Copper  tube, 

flow  rate,  however,  was  smaller  at  higher 

of  64- 

D,= 10  mm 

flow  rate;  hm  was  little  smaller  when  the 

95  kj  kg-1 

cold  carry  capacity  was  higher 

Song 

30  wt% 

Type  A  TBAB 

Stainless  steel 

0.35-6.8  ms-' 

Laminar  flow: 

hm  generally  increased  with  vs\  hm 

et  al. 

CHS,  xs=0- 

tube,  D,=14  mm, 

Num  =  1.231  X  10~5Re'  6606  Pr07073, 

increases  slightly  or  was  near  constant  in 

[27] 

18.4  wt% 

L= 2.8  m 

Turbulent  flow: 

the  turbulent  flow  region;  as  xs  increasing 

Num  =  5.254  x  10~4Re°9097  Pr1 1202 

to  a  certain  value,  hm  dropped  largely,  and 
the  flow  was  re-laminarized;  hm  increased 
with  xs  in  laminar  flow  region 

17.3  wt% 

Type  B  TBAB 

Laminar  flow: 

Variation  of  hm  was  similar  to  type  A  CHS, 

CHS,  xs=0- 

Num  =  1.798  x  10~3Re°6648  Pr13309, 

except  that  hm  in  turbulent  region 

26.2  wt% 

Turbulent  flow: 

Num  =  6.085  x  10~4Re°9,3°  Pr1  2061 

increased  rapidly  with  xs 

Ma  et  al. 

10  wt% 

Type  B  TBAB 

Copper  tube. 

0-2.5  ms-1  in 

Laminar  flow: 

h  was  enhanced  compared  to  water  at  the 

[28] 

CHS,  cps= 0- 

Dj—6,  14  mm, 

6  mm  pipe,  0- 

Nu(x)  =  0.2368(ReMR  Pr)° 5873  Ste-°  0795 

same  Reynolds  number;  h  decreased  at  the 

Ma  and 

Zhang 

[33] 


10  wt% 


Kumano  22.5  wt% 
et  al. 

[31] 

8.5  wt% 


20  vol% 


L=l.l  m 


1.2  m  s  1  in 
14  mm  tube 


A 


Transitional  and  turbulent  flow: 

Nu(x)  =  0.0219(ReMfi  Pr)07093Ste-°0498 

sq  \  0.2009  /A  \  -0.2441 


Type  B  TBAB 
CHS:  cps= 0- 
17.5  vol% 


Type  A  TBAB 
CHS,  0- 
20  wt% 

Type  B 
TBAB  CHS, 
0-20  wt  % 


Plate  heat 
exchanger,  10 
plates,  Dh= 4  mm, 
W =72  mm, 
H=310  mm 


Stainless  steel 
pipe,  D,  =  7.5, 
10.2  mm,  1=1  m 


2.5-13  L  min 


)Vu,oco,  =  0.338Re°,f7Pr'/3(i- 


AcoAH\ 

c^rJ 


1 000  <  Re  <  7500  Type  A  TBAB  CHS: 

Nu  =  3.33  Czft18x°42  (D//)° 
Type  B  TBAB  CHS: 

Nu  =  6.35  GzJi22x?  52  (D//)° 


Suziki  35  wt% 
et  al.  ammonium 

[157]  alum 


Copper  tube, 
D,  =  13  mm, 


12  wt% 
ammonium 
alum  hydrate  L=2  m 
slurry 


0.5-1.5  m  s 


melting  region  and  increased  at  the 
following  solution  region;  h  increased  with 
(ps  and  Reynolds  number,  while  cps  was 
more  important  for  laminar  flow  and 
Reynolds  number  dominated  the  heat 
transfer  in  turbulent  flow;  the  influence  of 
heating  power  was  distinct  in  the  case  of 
6  mm  tube 

hiocai  fluctuated  before  the  transition  from 
slurry  to  aqueous  solution,  and  increased 
sharply  when  turning  into  solution;  hoveraii 
basically  increased  with  the  inlet  water 
temperature,  however,  it  dropped  as  the 
increase  of  cps  until  an  un-fully  melting 
happened 

In  laminar  flow  region,  h  increased  with  xs, 
but  decreased  along  the  heating  tube, 
especially  at  the  tube  inlet,  and  the  effect 
of  Re  was  not  significant,  the  ratio  of 
measured  Nu  to  the  calculated  single¬ 
phase  Nu  was  about  1-2.5  and  increased 
with  xs;  in  turbulent  flow  region,  h  was 
almost  constant  with  xs,  but  dropped  at 
about  1 5  wt%,  then  increased  slightly  at 
high  xs,  h  changed  little  along  the  heating 
tube  and  increased  with  Re  distinctly,  the 
ratio  of  Nu  was  almost  1  with  all  xs 

With  the  same  Re,  Nu  of  slurry  was  higher 
than  that  of  solution  and  water;  Nu 
decreased  as  the  addition  of  surfactant 
caused  by  the  laminarization 


1.5-2. 5  times  heat  transfer  coefficient  was  achieved  in  compar¬ 
ison  with  smooth  tube.  The  same  conclusion  was  reached  by 
Alvarado  et  al.  [11]  about  MPCS.  They  compared  the  heat  transfer 
coefficient  of  MPCS  flowing  through  both  regular  and  enhanced 
8  mm  tubes  at  the  same  flow  velocity,  and  found  that  the  average 
heat  transfer  enhancement  was  about  30%  with  7  wt%  slurry. 
However,  the  enhancement  reduced  when  the  slurry  fraction 
increased  due  to  the  turbulence  damping. 


3.4.4.  Heat  transfer  characteristics  of  PCS  in  heat  exchangers 

In  addition  to  in  the  heated  tubes,  heat  transfer  performances 
of  PCS  in  the  heat  exchangers  have  also  been  investigated. 
Fig.  29(a)-(d)  shows  the  typical  heat  transfer  coefficients  of  ice 
slurry  in  the  plate  heat  exchangers  and  double-tube  heat  exchan¬ 
gers,  all  with  hot  water  in  the  other  side  of  heat  exchangers.  As 


the  increase  of  ice  fraction  from  5  to  22  wt%,  Bellas  et  al.  [135] 
obtained  nearly  constant  values  of  the  overall  heat  transfer 
coefficient  in  a  plate  heat  exchanger,  as  shown  in  Fig.  29(a). 
Moreover,  the  overall  heat  transfer  coefficient  reasonably 
increased  with  the  flow  rate.  The  authors  believed  that  the  overall 
heat  transfer  coefficient  was  the  function  of  ice  crystal  size, 
Reynolds  number  and  the  flow  passage  geometry.  Different  from 
the  results  reported  by  Bellas  et  al.  [135],  the  overall  heat  transfer 
coefficient  in  plate  heat  exchanger  measured  by  Norgaard  et  al. 
[136]  increased  with  the  inlet  ice  fraction,  as  shown  in  Fig.  29(b). 
As  shown  in  Fig.  29(c),  Lee  et  al.  [146]  found  that  the  heat  transfer 
coefficient  of  ice  slurry  flowing  through  double-tube  heat  exchan¬ 
ger  showed  weak  dependence  on  ice  fraction  when  the  flow  rate 
was  low,  and  there  was  a  sharp  increase  when  the  ice  fraction  was 
higher  than  about  1 0%.  However,  when  the  flow  rate  was  high,  the 
heat  transfer  coefficient  increased  gradually  as  the  ice  fraction 
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Fig.  23.  Ratio  of  Nusselt  number  of  ice  slurry  to  that  of  the  chilled  water  [112].  The  references  in  the  figure  are  those  in  [112]. 
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Fig.  24.  Local  heat  transfer  coefficients  and  Nusselt  numbers  of  ice  slurry  flowing 
through  heated  tubes  under  constant  heat  flux,  (a)  rectangular  channel, 
533  mm  x  100  mm  x  13  mm,  heated  at  one  side  by  hot  water  [112]  and 
(b)  generated  within  HYCOOL  fluid,  -35.5  “C,  21  mm  stainless  steel  pipe  [137], 


Fig.  25.  Local  heat  transfer  coefficients  of  MPCS  flowing  turbulently  through 
heated  tubes  under  constant  heat  flux,  (a)  Yamagishi  et  al.  [9],  10.1  mm  tube  and 
(b)  Wang  et  al.  [12],  4  mm  tube. 


increasing  from  0%  to  26%.  The  authors  also  pointed  that  the 
impact  of  hot  water  temperature  on  the  heat  transfer  coefficient 
could  be  neglected  in  spite  of  the  increase  of  heat  transfer  rate 
with  hot  water  temperature.  Bedecarrats  et  al.  [123]  presented 


the  variation  of  heat  transfer  coefficients  of  ice  slurry  flowing 
through  double-tube  heat  exchanger  with  smooth  and  corrugated 
inner  tubes,  and  the  typical  results  in  smooth  tube  is  shown  in 
Fig.  29(d).  It  was  noticed  that  the  variation  of  heat  transfer 
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Fig.  26.  Local  heat  transfer  coefficients  of  CHS  flowing  through  heated  tubes  under 
constant  heat  flux,  (a)  Kumano  et  al.  [31],  7.5  mm  tube  and  (b)  Ma  et  al.  [28], 
6  mm  tube. 

coefficient  was  similar  to  that  of  the  pressure  drop  shown  in 
Fig.  20(b).  At  lower  flow  rate  range,  the  flow  of  ice  slurry  with  0% 
ice  fraction  was  in  the  laminar  region,  therefore  the  increase  of  ice 
fraction  led  to  more  ice  particles  melting  near  the  heat  transfer 
wall  and  hence  the  heat  transfer  coefficient  increased;  at  higher 
flow  rate  range,  the  flow  of  ice  slurry  with  0%  ice  fraction  was  in 
the  turbulent  region,  and  the  turbulence  decreased  with  ice 
fraction  until  the  flow  changed  to  laminar  flow,  therefore  it  was 
noticed  that  the  heat  transfer  coefficient  increased  with  ice 
fraction  in  the  turbulent  region  and  followed  by  a  reduction  in 
the  transition  region,  and  then  a  more  rapid  increase  in  the 
laminar  region,  as  shown  in  the  figure.  The  authors  attributed 
the  rapid  increase  of  heat  transfer  coefficient  in  laminar  region  to 
the  thinner  ice-free  thermal  boundary  layer  which  was  caused  by 
the  enhanced  particle-wall  interactions  with  higher  ice  fraction. 
Higher  hot  water  temperature  was  found  to  result  in  higher  heat 
transfer  coefficient  both  in  laminar  and  turbulent  regions  except 
in  the  transitional  region.  Fernandez-Seara  et  al.  [118]  investi¬ 
gated  the  heat  transfer  performance  of  ice  slurry  in  a  terminal 
fan-coil  unit.  Results  showed  the  heat  transfer  was  enhanced  as 
using  ice  slurry  instead  of  chilled  water  by  the  factor  of  3.7-4.6, 
which  depended  on  the  fan  rotation  rate  as  well  as  the  ice 
fraction.  The  dominated  thermal  resistance  of  the  overall  heat 
transfer  was  on  the  air  side,  where  the  ratio  was  in  the  range  of 
80-88%.  However,  Ilian  and  Viedma  [158]  did  not  recommend  the 
direct  application  of  ice  slurry  in  fan-coil  unit,  since  the  frost 
formation,  the  partial  blockage  of  fan  coil  inlet  collector  and  the 
excessive  diameter  ratio  (dp/D,-)  displayed  the  disadvantages  of  ice 
slurry  compared  to  the  single-phase  flow. 

Ma  and  Zhang  [33]  investigated  the  heat  transfer  of  TBAB 
CHS  in  a  plate  heat  exchanger,  and  concluded  that  the  local  heat 


transfer  coefficient  was  highly  dependent  on  the  flow  velocity 
but  slightly  dependent  on  the  slurry  fraction,  and  /ilocal  would 
increased  drastically  after  the  slurry  turning  to  solution  due  to 
the  low  viscosity  and  high  Re  of  solution.  Meanwhile,  the 
overall  heat  transfer  coefficient  decreased  with  the  slurry 
fraction  both  at  a  certain  flow  velocity  and  inlet  hot  water 
temperature. 

3.4.5.  Heat  transfer  correlations 

Except  Knodel  et  al.  [132]  who  developed  the  heat  transfer 
correlation  of  ice  slurry  flowing  through  heated  tube  simply  by 
multiplying  a  factor  on  Petukhov  correlation  of  single-phase  fluid, 
other  researchers  considered  the  heat  transfer  coefficient  or 
Nusselt  number  of  ice  slurry  flowing  through  heated  tube  or  heat 
exchanger  as  the  complicated  function  of  Reynolds  number, 
Prandtl  number,  ice  fraction,  ratio  of  the  particle  diameter  to  the 
size  of  flow  passage  and  the  ratio  of  the  viscosity  of  bulk  fluid  to 
that  of  the  fluid  at  the  wall  and  so  on.  For  the  correlation  of  local 
Nusselt  number,  Graetz  number  {Gz=RePrDhjx)  is  widely  utilized 
[1 12,147-149],  The  index  of  Gz  in  the  correlation  is  positive,  which 
indicates  the  increment  of  Nusselt  number  with  Reynolds  number, 
Prandtl  number  as  well  as  the  entrance  distance.  The  index  of  ice 
fraction  is  also  positive  to  reflect  the  impact  of  ice  particles  on  heat 
transfer.  As  aforementioned,  the  influence  of  tube  diameter  on  heat 
transfer  lies  in  both  the  micro-convection  effect  of  ice  particles  and 
the  near-wall  thermal  layer,  therefore  the  index  of  diameter  ratio 
(dp/D,)  might  be  positive  and  negative  [133,145,148,150]  under 
different  situations.  Niezgoda-Zelasko  [145]  and  Niezgoda-Zelasko 
and  Zelasko  [148]  recommended  a  term  of  ((AH/CpAr)(Axs/100)) 
to  reflect  the  influence  of  ice  melting  on  the  heat  transfer  when 
developing  the  heat  transfer  correlation,  which  showed  good 
agreement  with  experimental  results. 

Similar  to  ((AH/CpAT)(Axs/100)),  Charunyakorn  et  al.  [159] 
proposed  the  utilization  of  Stefan  number,  (Ste=  (CpqwDJ 
(2As))/(xsAH)),  to  reflect  the  influences  of  latent  heat  and  heating 
power  on  the  heat  transfer  in  developing  the  heat  transfer 
correlation  of  MPCS  flowing  through  heated  tube,  and  Goel  et  al. 
[151]  validated  the  importance  of  Ste  on  the  heat  transfer  of  MPCS 
experimentally,  therefore  Wang  et  al.  [13]  used  this  parameter  to 
develop  the  heat  transfer  correlation  of  MPCS  in  heated  tube  and 
the  relative  error  was  shown  to  be  within  + 10%.  Ma  et  al.  [28] 
also  utilized  this  term  for  TBAB  CHS  flowing  through  heated  tube 
and  the  similar  term  of  (1  +AxsAH/(CpAT))  for  TBAB  CHS  flowing 
through  plate  heat  exchanger  [33].  Good  results  were  achieved 
using  these  latent  heat  terms. 

4.  Application  of  PCS  to  refrigeration  and  air  conditioning 
and  performance  evaluation 

4.3.  Ice  slurry 

Because  of  the  numerous  advantages,  ice  slurry  is  undoubtedly 
promising  and  can  be  used  in  various  applications,  including  air 
conditioning,  refrigeration,  process  cooling  and  so  on.  Many  ice 
slurry  systems  were  installed  in  Japan,  which  mainly  concen¬ 
trated  in  air  conditioning  system  with  thermal  energy  storage, 
while  the  applications  in  Europe  and  USA  cover  all  the  fields.  The 
typical  utilization  of  ice  slurry  can  be  divided  into  the  following 
ways:  direct  contact  cooling  application,  typically  represented  by 
fish  cooling  system;  direct  application,  ice  slurry  flows  directly  to 
the  load  side  and  exchanges  the  thermal  energy  with  the  high 
temperature  air  or  water;  direct  carrier  liquid  application,  only 
liquid  phase  in  the  slurry  is  pumped  to  the  load  side  to  avoid  the 
blockage  in  the  narrow  flow  passages;  indirect  application,  ice 
slurry  is  only  employed  as  the  thermal  storage  material  while  the 


P.  Zhang,  Z.W.  Ma  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  5021-5058 


5051 


1.02 

1.00 

0.98 

0.96 

0.94 

0.92 

0.90 

0.88 

0.86 

0.84 


I 

♦ 


V* 

V* 


%  ui 

9m  . 


Re  =  55000 


_L_ 


0  2  4  6  8  10 

ICE  FRACTION,  %  by  mass 


12 


c  d 


Fig.  27.  Mean  heat  transfer  coefficients  and  Nusselt  numbers  of  ice  slurry  flowing  through  heated  tubes  under  constant  heat  flux,  (a)  24  mm  stainless  steel  tube, 
q=15  kW  m~2  [132];  (b)  16  mm  cooper  tube,  w  is  flow  rate  (m  s_1)  and  q  is  heat  flux  (W  m~2)  [145];  (c)  21  mm  stainless  steel  tube  [149]  and  (d)  smooth  stainless  steel 
tube,  q>  is  the  inlet  ice  fraction  [133J. 


cold  energy  is  first  exchanged  and  delivered  to  the  load  side  by 
another  coolant,  such  as  chilled  water. 

A  typical  application  of  direct  ice  slurry  application  system  in 
large-scale  building  is  the  CAPCOM  building  in  Osaka,  Japan 
[3,160],  The  building  is  20  stores  and  totally  has  16,784  m2  floor 
area,  where  the  2nd  to  1 6th  floor  are  cooled  by  ice  slurry.  The  ice 
slurry  is  generated  and  stored  at  the  17th  floor,  as  shown  in 
Fig.  30.  There  are  two  ice  slurry  units  driven  by  gas  driven 
absorption  refrigeration  system  with  272  kW  cooling  capacity 
each.  The  ice  slurry  is  generated  by  a  shell-and-tube  ice  generator, 
and  ice  is  formed  inside  the  tube  and  is  then  stored  in  96  m3 
storage  tank.  Each  unit  has  three  alternating  operated  heat 
exchangers,  and  two  for  freezing  and  one  for  defrosting.  8  mixers, 
which  are  totally  55  kW,  are  installed  in  the  tank  to  agitate  the 
mixture  of  ice  and  water,  then  40  wt%  ice  slurry  at  the  tempera¬ 
ture  range  of  -2.4  to  -2.6  °C  is  obtained.  The  ice  slurry  is 
pumped  from  the  tank  first  through  an  ice  slurry  separator  to 
adjust  the  slurry  fraction  to  20  wt%,  and  then  reaches  the  air 
handing  units  on  the  floors  of  2-16.  Due  to  the  low  temperature 
and  large  cooling  capacity  of  ice  slurry,  the  air  distribution 
temperature  in  the  handing  unit  is  designed  to  12  °C  which  is 
lower  than  the  normally-used  temperature  at  1 5  °C,  therefore  the 
air  flow  rate  is  reduced  by  about  22%  which  leads  to  capital  and 
operation  cost  reduction.  Compared  to  the  predicted  performance 


of  an  equivalent  conventional  ice  thermal  storage  system  which 
has  higher  air  distribution  temperature,  the  total  energy  con¬ 
sumption  of  ice  slurry  system  during  the  summer  of  1996 
achieved  4%  reduction  resulted  from  the  energy  saving  in 
hydraulic  and  air  distribution. 

An  indirect  ice  slurry  system  was  built  for  a  190,000-square- 
foot  basketball  arena  and  athletic  complex  cooling  usage  at 
Virginia  Commonwealth  University  in  USA  [161],  which  has  a 
seating  capacity  of  7500  people.  The  system  schematic  is  shown 
in  Fig.  31.  The  base  aqueous  solution  used  for  ice  slurry  genera¬ 
tion  is  7%  glycol  solution,  of  which  the  freezing  temperature  is 
-2.2  °C.  As  shown  in  the  figure,  the  solution  is  pumped  to  the 
ice  slurry  generator  which  has  the  capacity  of  380  t  (about 
1336  kW),  and  -8.3  °C  primary  refrigerant  is  utilized  for  ice 
generation.  The  ice  slurry  is  generated  and  stored  at  -3.3  °C. 
However,  the  ice  slurry  is  not  pumped  to  the  load  side  directly, 
only  the  chilled  solution  at  1.1  °C  is  extracted  from  the  tank  to 
the  building  heat  exchanger  to  cool  water,  which  is  located 
about  100  yards  from  the  ice  slurry  system.  After  heat  exchan¬ 
ging,  warm  solution  at  11.6  °C  returns.  Thereafter,  2.2  °C  water  is 
supplied  to  the  air  handing  unit  to  satisfy  the  proper  fresh  air  in 
the  7500  seat  arena.  Compared  to  the  conventional  5-7  DC 
coolant,  the  supply  of  2.2  °C  water  causes  the  scale  reduction 
of  many  devices,  such  as  air  coil,  fan  motor,  water  pump,  duct 
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Fig.  28.  Mean  heat  transfer  coefficients  of  type  ATBAB  CHS  flowing  through  heated  tubes  under  constant  heat  flux,  (a)  14  mm  tube.  Song  et  al.  [27]  and  (bl),  (b2)  7.5  mm 
tube,  Kumano  et  al.  [31]. 


and  pipe  etc.  The  annual  operation  cost  is  reduced  by  approxi¬ 
mately  $75,000. 

Sunwell  Technologies  Inc.  [162]  has  supplied  a  number  of  direct 
contact  ice  slurry  systems  for  fish  processing  and  preservation  as 
well  as  for  supermarket  display  usage.  The  named  Deepchill 
Variable-State  Ice  System  can  provide  0-100  wt%  ice  slurry  with 
0.25-0.5  mm  round  crystals  at  the  temperature  range  of  -  30  to 
0  °C,  which  is  easily  pumpable.  All  components  that  employed  in 
the  ice  slurry  generation,  storage  and  delivery  are  constructed  by 
medical  grade  stainless  steel  or  food  grade  plastic  composites, 
while  the  slurry  is  always  sealed  and  never  exposed  to  outside 
contaminants.  Hence,  the  slurry  is  considered  as  clean  and  reliable 
enough  for  direct  contact  with  food.  25-40  wt%  ice  slurry  is  usually 
utilized  and  can  cool  the  target  goods  to  -  2  to  0  °C  in  several 
minutes. 

Ice  slurry  system  constructed  in  a  large  European  institutional 
kitchen  in  France  [4,163]  provides  a  complete  cooling  chain  in 
this  kitchen,  including  refrigerated  room,  refrigeration  system, 
refrigerators,  recooling  vessels,  the  food  service  conveyor  belts, 
plate  dispensers  and  cooling  tray  transport  trolleys.  This  system 
employs  16  scraped  surface  ice  slurry  generator  of  185  kW  total 
cooling  capacity  and  3  cylindrical  storage  tank  of  22  m3  to 
generate  and  store  ethanol/water  ice  slurry.  It  was  found  that 
the  usage  of  ice  slurry  could  provide  a  constant  low  temperature 
for  food  processing  and  satisfy  all  the  cooling  requirements  in  the 
large  kitchen 

4.2.  MPCS 

An  MPCS  thermal  energy  storage  system  was  constructed  for 
cooling  in  Narita  Airport  in  Japan  [18]  without  replacing  the 
existed  2500  RT  (about  8792.5  kW)  electrical  turbo  chiller.  This 
project  aimed  to  improve  the  decrease  in  capacity  and  perfor¬ 
mance  caused  by  refrigerant  substituting  for  environmental  con¬ 
servation.  Non-toxic  and  phase  change  temperature  controllable 


material,  n-Paraffin  waxes,  was  selected  as  the  PCM.  The  slurry  is 
the  mixture  of  water  and  2  pm  microcapsules,  and  the  melting 
temperature  and  the  corresponding  latent  heat  of  which  are  8  °C 
and  75.9  kjkg-1,  respectively.  The  storage  tank  is  24.7  m  in 
height  and  7.4  m  in  diameter.  Cold  energy  is  stored  during  night 
time.  MPCS  at  11.5  °C  is  pumped  by  ordinary  centrifugal  pumps 
from  the  top  of  the  storage  tank,  and  is  cooled  in  a  plate-fin  type 
heat  exchanger  by  3.5  °C  coolant,  as  shown  in  Fig.  32(a).  Phase 
change  process  occurs,  the  slurry  is  then  cooled  to  4  °C  and  stored 
at  the  bottom  of  the  tank.  The  stored  cold  energy  is  released 
during  day  time  as  shown  in  Fig.  32(b).  By  the  heat  exchange  with 
cold  MPCS,  the  coolant  of  5.1  °C  is  obtained  and  pumped  to  the 
load  side.  COP  of  the  chiller  and  the  system  during  the  cold 
storage  process  is  about  5.4  and  3.2,  respectively,  which  is  lower 
than  the  COP  of  5.9  and  4.4  when  there  is  no  energy  storage. 
Compared  to  the  old  storage  system  using  external  melting  ice, 
though  the  thermal  storage  density  of  the  used  MPCS  system  is 
about  60%  smaller,  the  operation  cost  is  32%  lower. 

Griffiths  and  Eames  [16]  investigated  the  performance  of  MPCS 
as  the  heat  transfer  medium  in  a  test  chamber  containing  a 
chilled  ceiling  panel,  and  the  schematic  diagram  of  which  is 
shown  in  Fig.  33.  The  authors  employed  the  MPCS  with  a  melting 
temperature  of  18  °C  and  2-8  pm  microcapsules  produced  by 
BASF,  and  the  phase  change  temperature  is  close  to  the  working 
temperature  of  chilled  ceiling  at  16-18  °C.  The  authors  compared 
the  performance  of  chilled  water  and  MPCS  with  40%  concentra¬ 
tion.  To  maintain  the  test  chamber  at  a  constant  temperature  of 
19  °C,  flow  rate  of  0.7  L  s_1  was  required  when  water  was  utilized 
as  the  heat  transfer  medium,  while  only  0.25  Ls~]  was  required 
for  40%  MPCS  to  maintain  the  space  temperature  at  20-21  °C, 
since  the  average  specific  heat  of  the  slurry  in  the  temperature  of 
16-17.5  °C  was  as  high  as  9  k]  kg-1  K_1,  with  sometimes  a  value 
of  11-12  kj  kg-1  K-1.  Through  the  continuous  operation  for  over 
four  months,  the  MPCS  with  40%  concentration  was  validated  as 
heat  transfer  medium  on  chilled  ceiling  panel  instead  of  water, 
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Fig.  29.  Overall  heat  transfer  coefficients  of  ice  slurry  through  heat  exchangers,  (a)  plate  heat  exchanger,  24  plates,  Dh= 4  mm,  IV—  1 1 2  mm,  H=310  mm  [135];  (b)  plate 
heat  exchanger,  10  plates,  Dh= 4  mm,  IV— 70  mm,  H= 300  mm  [136];  (c)  double-tube  heat  exchanger,  inner  tube  D,= 13.84  mm,  outer  tube  Dt= 25  mm,  l=l.5m  [146]  and 
(d)  double-tube  heat  exchanger,  inner  tube  D,= 22  mm,  outer  tube  Dt= 33  mm,  1=2  m  [123], 


and  there  was  neither  deposition  of  microcapsules  in  pipe  nor  the 
degradation  of  the  slurry,  and  the  system  pumping  power  was 
reduced  due  to  small  flow  rate.  However,  the  authors  noticed 
that  the  control  procedure  for  cooling  the  slurry  should  be 
optimized. 

4.3.  TBAB  CHS 

Takao  et  al.  [21  ]  introduced  the  operation  results  of  the  TBAB 
CHS  air  conditioning  system  in  the  office  building  of  NKK  Keihin 
Works,  which  was  constructed  by  the  JFE  Engineering  Corpora¬ 
tion.  Total  area  of  the  two  floors  building  is  about  4000  m2,  where 
1700  m2  is  the  cooling  object  of  TBAB  CHS.  As  shown  in  the 
system  sketch  in  Fig.  34,  an  absorption  refrigerator  is  utilized.  The 
cold  carry  capacity  of  the  employed  TBAB  CHS  is  approximate 
2  times  of  that  of  chilled  water,  which  is  about  59  MJ  m  3, 
therefore  the  coolant  flow  rate  was  reduced  by  50%.  The  heat 
exchange  between  TBAB  CHS  and  air  inside  the  room  was  under¬ 
taken  by  plate  fin  heat  exchanger,  and  the  diameter  of  the  cooper 
tube  for  TBAB  CHS  was  reduced  to  10  mm  because  of  the  low  flow 
rate.  The  overall  heat  transfer  coefficient  in  the  heat  exchanger 
was  about  250-300  W  m-2  K  1  as  the  CHS  flow  rate  varying  from 
1  to  10  kg  min-1,  which  was  higher  than  that  of  water  when  the 


flow  rate  was  smaller  than  about  6.0  kg  min-1.  However,  the 
overall  heat  transfer  coefficient  was  lower  at  higher  flow  rate 
range.  Compared  to  the  utilization  of  chilled  water,  though  there 
was  a  20%  increase  of  the  pressure  drop  of  TBAB  CHS,  the 
operation  pumping  power  measured  in  the  representative  day 
of  August  of  2000  was  reduced  by  47%  due  to  the  lower  flow  rate, 
while  the  total  energy  consumption  had  a  reduction  by  21%. 

Ogi  [164]  reported  the  thermal  energy  storage  air  conditioning 
system  using  TBAB  CHS  in  Kawasaki  Underground  Azalea  which 
has  a  total  area  of  about  56,000  m2.  The  project  was  also 
conducted  by  JFE  Engineering  Corporation  and  was  finished  in 
2008.  The  system  schematic  diagram  is  shown  in  Fig.  35.  The 
system  was  modified  based  on  an  existing  air  conditioning 
system,  which  consisted  an  air  handling  unit  (AHU)  system  using 
3  heat  pump  of  550  Rt  (about  1934  kW)  each,  and  a  fan  coil  unit 
(FCU)  system  with  a  centrifugal  chiller  of  500  Rt  (about  1758  kW). 
Thermal  energy  storage  by  water  was  employed  by  the  FCU 
system  using  a  storage  tank  of  700  m3.  When  TBAB  CHS  was 
introduced  to  the  FCU  system  instead  of  water,  the  centrifugal 
chiller  was  substituted  by  a  screw  chiller  with  500  Rt  cooling 
capacity  and  the  COP  was  improved  from  4.19  to  4.44.  Since  the 
cold  carry  capacity  of  TBAB  CHS  was  about  2  times  as  large  as  that 
of  the  chilled  water,  the  stored  cold  energy  was  improved  from 
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*  Values  in  (  )  indicate  those  at  the  second  and  final  stages. 
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Fig.  31.  Schematic  diagram  of  indirect  ice  slurry  application  in  university  arena 
[161]. 


840  Rth  (about  2954  kWh)  by  water  to  1650  Rth  (about 
5803  kWh)  by  TBAB  CHS  without  changing  the  storage  tank. 
Compared  to  the  original  water  storage  system,  the  new  TBAB 
CHS  storage  system  led  to  a  40.17%  reduction  of  electricity 
consumption  and  a  40.28%  reduction  of  C02  emission. 

Moreover,  all  the  other  TBAB  CHS  projects  conducted  by  JFE 
Engineering  Corporation,  including  the  systems  of  Takenaka 
Corporation  (2004),  Tsurumi  office  building  of  JFE  Engineering 
Corporation  (2005),  Kuraray  Kurashiki  Office  Production  Technol¬ 
ogy  Center  (2007),  NTN  Okayama  Plants  (2009),  achieved  more 
than  20%  energy  saving  [165,166]. 


Fig.  32.  Operation  strategies  of  the  MPCS  system  during  (a)  night  time  and  (b)  day 
time  [18]. 


5.  Conclusions 

In  this  paper,  the  application  of  phase  change  material  slurries, 
i.e.,  ice  slurry,  MPCS  as  well  as  CHS  in  the  secondary  loop 
refrigeration  and  air  conditioning  system  have  been  reviewed, 
and  the  utilizations  of  PCS  are  expected  to  facilitate  the  energy 
saving  and  environmental  conservation.  The  paper  mainly  focuses 
on  the  generation  method  and  storage  strategies,  the  flow  and 
heat  transfer  characteristics  and  the  applications  of  these  slurries 
in  real  projects.  There  are  many  common  generation  methods  of 
ice  slurry  in  which  the  necessary  technologies  are  used  to  prevent 
the  generated  solid  crystals  from  adhering  to  the  heat  transfer 
surface  to  deteriorate  the  heat  transfer,  e.g.,  using  scraper  or  rod 
or  fluidized  solid  particles  to  detach  the  adherence,  using  the 
supercooling  method  to  avoid  the  crystals  forming  on  heat 
transfer  surface,  using  direct  contact  heat  exchange  to  improve 
the  heat  transfer,  and  so  on.  The  producing  process  for  MPCS  is 
relatively  complicated  and  higher  cost,  which  may  limit  its 
application.  The  homogeneous  storage  of  ice  slurry  and  CHS 
benefits  the  extraction  but  has  limitation  on  storage  density, 
while  the  heterogeneous  storage  can  achieve  high  storage  density 
but  requires  additional  power  for  extraction  and  the  crystals  size 
may  become  larger.  Though  PCS  has  been  widely  investigated  and 
utilized,  the  fundamental  investigations  on  flow  and  heat  transfer 
characteristics  are  still  divergent  because  of  the  complex  influ¬ 
ences  by  many  factors,  therefore  further  studies  are  still  required. 
Moreover,  the  application  of  these  slurries  show  quite  positive 
significance  on  energy  saving  and  environment  conservation. 
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Fig.  33.  Schematic  diagram  of  chilled  ceiling  in  the  testing  chamber  [16]. 
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Fig.  34.  Schematic  diagram  of  air  conditioning  system  using  TBAB  CHS  in  the  office  building  of  NKK  Keihin  Works  [21]. 


Fig.  35.  Schematic  diagram  of  air  conditioning  system  using  TBAB  CHS  in  Kawasaki  Underground  Azalea  [164]. 
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However,  MPCS  and  CHS  are  still  not  thoroughly  investigated  and 
applied,  which  suggests  that  more  fundamental  researches  are 
still  necessary  for  MPCS  and  CHS. 
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